UNCLASSIFIED 


AD  NUMBER 


AD875549 


NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors;  Critical 
Technology;  APR  1970.  Other  requests  shall 
be  referred  to  School  of  Aerospace 
Medicine,  Brooks  AFB,  TX. 


AUTHORITY 


SAM  USAF  ltr ,  10  Sep  1976 


THIS  PAGE  IS  UNCLASSIFIED 


THIS  REPORT  HAS  BEEN  DELIMITED 
AND  CLEARED  FOR  PUBLIC  RELEASE 
UNDER  DOS  DIRECTIVE  5200,20  AND 

NO  RESTRICTIONS  ARE  IMPOSED  UPON 
ITS  USE  AND  DISCLOSURE, 

DISTRIBUTION  STATEMENT  A 

APPROVED  FOR  PUBLIC  REMEASE; 
DISTRIBUTION  UNLIMITED, 


fip  $*7  S' 


Final  Report  -  Pert  / 


ATMOSPHERIC  EFFECTS  UPON  LASER  EYE  SAFETY 


By:  WARREN  B.  JOHNSON  WILLIAM  E.  EVANS  EDWARD  E.  UTHE 


Prepared  for: 

UNITED  STATES  AIR  FORCE 
SCHOOL  OF  AEROSPACE  MEDICINE 
AEROSPACE  MEDICAL  DIVISION 
BROOKS  AIR  FORCE  BASE,  TEXAS 


D  D  C 

irainnE^ 

OCT  *0  1970 


STANFORD  RESEARCH  INSTITUTE 

ICanlo  Park,  California  Mon  •  u.s.A. 


Reproduced  by  the 

CLEARINGHOUSE 
for  Federel  Scientific  &  Technical 
Information  Springfield  Va.  22151 


Final  Report  -  Part  I 


April  1970 


ATMOSPHERIC  EFFECTS  UPON  LASER  EYE  SAFETY 


a 


tn 

§ 


x 

E> 


% 


H 

M 


1  "  |  By:  WARREN  B.  JOHNSON  WILLIAM  E.  EVANS 


n  o  | 

H  rt  ! 

a  .  k  * 


CS  V 1 


Prepared  for: 


UNITED  STATES  AIR  FORCE 
SCHOOL  OF  AEROSPACE  MEDICINE 
,  AEROSPACE  MEDICAL  DIVISION 
^JBROOKS  AIR  FORCE  BASE.  TEXAS 

CONTRACT  F41 609-69-C-0001 


♦> 

o  C  u 

9  0  0 

JSh 

L.  e) 
o  ®  > 
►  o 
o  h 
M IQ. 

o  a 

ecu 

-*~i  GO 

.o  h 

3  ii  O 
cl  t.  -< 
o  u 
m  v.  o, 


SRI  Project  7472 


Approved  by: 

R.  T.  H.  COLLIS.  Director 
Aaropbytk*  Laboratory 


o 

*>  *» 
Q 

§2 
o  •> 

O  ■*« 
^0  8 

01 

w  r: 
•<  ® 
JC  t 

»-•  w 


RAY  L.  LEAOABRAND.  executive  Director 
electronic*  and  Radio  Science*  DMdon 


EDWARD  E.  UTHE 


Copy  No. 


ABSTRACT 


The  results  of  cwo  laser  scintillation' experiments  in  the  atmosphere 
are  presented,  as  is  a  revised  eye-hazard  evaluation  procedure  incorpo¬ 
rating  these  results.  The  objective  of  the  research  is  to  determine  the 
influence  of  atmospheric  thermal  turbulence  on  laser  eye  safety. 

The  experimental  work  on  a  3.5-km  near-ground  horizontal  propagation 
path  was  designed  to  determine  the  effect  of  range  upon  laser  scintilla¬ 
tion  magnitude,  represented  by  the  log  intensity  standard  deviation  (c). 
It  was  found  that  for  very  short  ranges  and  low  thermal  turbulence,  the 

measured  values  (C  )  agree  well  with  theory  (a  ).  However,  for  longer 
m  t 

ranges  and  higher  thermal  turbulence,  the  values  of  a  reached  a  maximum 

m 

value  of  approximately  1.25  at  a  =  2.50,  and  then  decreased  beyond, 

2 

fitting  the  empirical  function  a  =  a  /(l  +  0.16  a  )  fairly  well.  This 

m  t  t 

"supersaturation  effect"  is  quite  significant  in  reducing  eye-damage 
probabilities  for  typical  ranges  and  thermal  turbulence  conditions. 

The  second  experiment  was  an  attempt  to  confirm  that  scintillation 
saturation  occurs  on  a  slant  path  between  a  ground-based  laser  and  a 
receiver  on  an  orbiting  aircraft.  Some  evidence  of  saturation  was  found, 
but  the  effect  was  not  clearly  established  because  of  the  limited  ranges 
necessitated  by  the  laser  used. 

The  eye-hazard  evaluation  procedure  is  quite  generally  applicable, 
and  requires  only  readily  available  input  data.  The  method  covers 
propagation  on  a  slant  path  with  either  a  ground-based  or  airborne 
laser,  as  well  as  near-ground  horizontal  propagation.  The  airborne- 
laser  slant-path  configuration  is  shown  to  give  the  lowest  eye-damage 
probability,  other  factors  being  the  same. 
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I  INTRODUCTION 


] 

1 

! 

1 

i 

I 

The  rapid  proliferation  of  laser  systems  for  civilian  and  military  I 

applications  urgently  requires  the  development  of  suitable  guidelines  j 

I 

1 

to  ensure  safety  in  their  use.  In  the  absence  of  appropriate  and  ] 

realistic  guidelines,  excessively  conservative  "safety  factors"  may  be  j 

proposed  that  could  unduly  restrict  the  use  of  laser  systems.  I 

The  U.S.  Air  Force  is  actively  conducting  and  sponsoring  research  < 

studies  (1)  to  determine  threshold  levels  for  laser  eye  damage 

*  t 

(Vassiliadis  et  al.,  1969),  and  (2)  to  determine  the  effects  of  atmo¬ 
spheric  thermal  turbulence  on  eye-damage  probabilities.  The  latter  is 
the  objective  of  our  present  work  and  the  subject  of  this  report.  A 
previous  report  (Johnson  et  al.,  1968)  described  the  results  of  the 
first  year's  research. 

The  effect  of  atmospheric  temperature  inhomoger.eities  in  increasing 
eye-damage  probabilities  is  exemplified  in  Figure  1,  which  illustrates 
the  resulting  characteristic  breakup  of  a  propagating  laser  beam.  To  a 
viewer,  this  appears  as  scintillation  (fluctuations  in  intensity).  The 
"hot  spots,"  or  areas  of  localized  beam  intensification,  typically  have 
intensities  that  are  tens  or  hundreds  of  times  the  average  beam  intensity. 

Thus,  as  the  state  of  the  atmosphere  changes  from  a  thermally  "quiet" 
situation  (no  scintillation)  to  ore  of  increasing  thermal  turbulence  and 
scintillation,  the  probability  of  laser  eye  damage  increases  drastically. 

During  the  first  year,  substantial  progress  was  made  toward  clari- 
lying  and  parameterizing  the  effects  of  atmospheric  thermal  turbulence 


« 

References  are  listed  at  the  end  of  the  report. 
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FIGURE  t  BEAM  PATTERN  OF  RUBY  LASER  AT  11.6  km 
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upon  laser  propagation  and  upon  the  associated  eye-damage  probabilities. 

It  was  determined  experimentally  that  the  amplitude  probability  distri¬ 
bution  of  the  laser  intensity  measured  by  a  5-mm-aperture  detector  is 
well  described  by  a  log-normal  distribution,  in  accordance  with  Tatarski’s 
(1961)  theory.  The  results  demonstrated  that,  although  the  observed 
probability  distributions  departed  somewhat  from  log  normal  at  the  ex¬ 
treme  upper  end  of  the  distribution  (possibly  because  of  the  finite 
laser  power),  the  assumption  of  a  log-normal  distribution  gave  safe 

(conservative)  e>e- damage  probability  estimates.  These  analyses  were 

-7 

carried  out  to  probabilities  approaching  10  ,  almost  three  orders  of 

magnitude  lower  than  previously  available. 

Using  the  log-normal  model,  a  preliminary  set  of  laser  safety  guide¬ 
lines  for  atmospheric  effects  were  developed.  These  tentative  guidelines 
are  based  upon  readily  available,  conventional  weather  observations. 
However,  a  number  of  approximations  had  to  be  made  in  developing  the 
guidelines,  because  of  the  basic  lack  of  information.  During  the  past 
year,  additional  experiments  were  carried  out  in  order  to  clarify  certain 
unresolved  propagation  effects  and  to  permit  the  guidelines  to  be  extended 
to  a  broader  range  of  usage  of  laser  systems  in  terms  of  path  geometry, 
nature  of  the  platform,  and  turbulence  conditions.  The  results  of  this 
research  and  the  incorporation  of  this  information  into  the  laser  eye- 
safety  guidelines  are  described  in  the  remainder  of  this  report.  It 
should  be  emphasized  that  these  safety  guidelines  should  be  considered 
tentative  and  used  with  caution  until  the  new  experimental  results  that 
they  reflect  can  be  confirmed  for  a  variety  of  laser-beam  and  propagation- 
path  configurations. 
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II  TECHNICAL  BACKGROUND,  SUMMARY  OF  RESULTS,  AND  RECOMMENDATIONS 


A.  Background  and  Summary  of  Results 

Classical  optical  propagation  theory  (Tatarski,  1961)  predicts  that 

the  scintillation  magnitude  (a^)  is  directly  proportional  to  thermal 

turbulence  intensity  (C  ),  approximately  directly  proportional  to 

n 

range  (R),  and  approximately  inversely  proportional  to  the  square  root 
of  the  optical  wavelength;  the  exact  form  is 

11/12  ,  v7/12 

O  -  K  R  (2ttA)  c  .  (1) 

t  n 


Here  o  is  the  theoretical  log  signal  intensity  standard  deviation,  and 
K  is  a  dimensionless  constant.  This  theory,  which  is  based  on  the  so- 
called  Rytov  approximation,  fits  observations  well  for  low-thermal- 
turbulence  conditions  and  short  ranges. 

Application  of  the  theory  to  longer  ranges  and  higher  turbulence 
levels  has  not  been  intuitively  appealing,  however,  since  for  these  con¬ 
ditions  the  predicted  scintillation  magnitude  (and  hence  the  eye-damage 
probability)  increases  indefinitely  without  bound.  But  no  better  theory 
has  been  available,  and  indeed,  recent  advances  have  come  principally 
through  experimentation. 


The  first  breakthrough  came  from  experiments  in  Russia  by  Gracheva 
and  Gurvich  (1965)  and  by  Gracheva  (1967),  which  showed  a  saturation, 
or  levelling  off,  of  the  scintillation  magnitude  with  increasing  range 
beyond  about  1  km.  (The  longest  range  at  which  measurements  are  reported 
in  the  Russian  work  is  1.75  km.)  Subsequently,  Tatarski  (1966)  and 
DeWolf  (1968)  developed  variations  of  the  propagation  theory  that  con¬ 
form  to  some  degree  to  the  Russian  observations. 
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In  November  1968,  at  an  excellent  site  near  Woodland,  California, 
we  carried  out  a  horizontal  propagation  experiment  designed  primarily  to 
check  and  extend  the  Russian  results  regarding  the  range  dependence  of 
the  saturation  magnitude.  Our  results,  which  cover  a  broad  range  of 
thermal-turbulence  conditions  and  ranges  out  to  3.5  km,  were  surprising 
and  are  believed  to  be  extremely  significant  in  the  context  of  the  eye- 
safety  problem.  We  found  that,  for  our  experimental  conditions,  the 
scintillation  magnitude  not  only  saturated  with  increasing  range  and 
thermal  turbulence  level  (as  represented  by  a^),  but  actually  reached  a 
peak  and  decreased  thereaf ter — an  effect  that  we  have  called  "super¬ 
saturation."  An  analytical  expression  that  fits  our  measured  data  well 
is 

am  -  - “ - 2  >  <2' 

1  +0.16  a 

where  a  is  the  measured  log  intensity  standard  deviation, 
m 

Subsequent  to  the  horizontal  propagation  experiment,  we  learned 
that  Deitz  and  Wright  (1969)  had  reported  some  indications  of  such  an 
effect  in  their  data,  which  extended  in  range  to  1.5  km.  In  addition, 
Mevers  et  al.  (1969)  reported  in  April  1969  that  their  data,  which 
extended  to  5  km,  definitely  shows  a  supersaturation  effect;  Ochs  (1969) 
also  hints  of  such  an  effect  in  his  results.  Thus,  three  or  four  indepen¬ 
dent  experiments,  in  which  the  investigators  were  unaware  of  each  other's 
results,  have  produced  results  that  indicate  the  existence  of  the  super¬ 
saturation  phenomenon.  '  Furthermore,  a  close  examination  of  the  Russian 
data  reveals  that  they,  too,  reveal  such  a  tendency  within  the  limited 
range  covered. 

Other  significant  results  obtained  from  the  horizontal  propagation 
experiment  include  the  following: 
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(1)  The  frequency  spectrum  of  the  scintillation  shifts  to 
higher  frequencies  with  longer  range  and  increasing 
supersaturation,  and  tends  to  become  bimodal. 

(2)  The  intensity  probability  distribution  remains  approxi¬ 
mately  log-normal  for  both  subsaturation  and  super¬ 
saturation  conditions,  for  a  detector  centered  in  the 
laser  beam. 

(3)  Systematic  deviations  from  log-normal  occur  when  the 
detector  is  not  centered  in  the  laser  beam. 

(4)  Values  of  C  derived  from  measurements  with  a  fast- 

n 

response  differential  thermometer  system  show  excellent 
correlation  with  scintillation  magnitude  for  subsatura¬ 
tion  conditions,  in  accordance  with  Eq .  (1).  The 
constant  K  was  determined  to  be  1.06. 

A  slant-path  experiment  using  an  orbiting  aircraft  was  carried 
out  in  April  1969  at  Brooks  AFB.  The  quantity  of  data  obtained  (28  good 
runs  in  two  days)  was  less  than  desired  because  of  the  limited  availa¬ 
bility  of  the  aircraft.  Some  indications  of  saturation  of  scintillation 
on  the  slant  path  were  found  in  the  limited  data  sample  obtained.  The 
maximum  range  used  was  approximately  4000  m,  which  was  not  far  enough 
to  produce  supersaturation.  Saturation  effects  were  observed  on  the 
horizontal  control  range  at  a  distance  of  350  m. 

The  pseudo-wind  effects  caused  by  the  aircraft  and  laser-beam  move¬ 
ment  were  readily  apparent  in  the  higher  frequencies  of  scintillation 
on  the  slant  path,  typically  greater  than  1  kHz. 
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B.  Recommendations  for  Additional  Research 


Although  available  information  shows  that  the  supersaturation  effect 

definitely  exists,  the  data  from  the  experiments  by  various  investigators 

do  not  agree  in  terms  of  the  maximum  scintillation  magnitude  and  the  rate 

of  decrease  of  scintillation  with  increasing  range  and  thermal  turbulence 

* 

level.  The  maximum  scintillation  magnitude  (O  )  effectively  sets  an 

m 

upper  limit  for  the  eye-damage  hazard  (probability).  Thus,  it  is  partic- 

★ 

ularly  important  for  the  eye-safety  problem  that  a  be  accurately 

m 

determined,  and  that  we  determine  the  nature  of  any  dependence  upon 
propagation-path  geometry  (e.g.,  elevation),  meteorological  conditions, 
or  laser  characteristics.  Additional  measurements  on  the  slant  path 
and  also  at  night  under  stable  conditions  are  needed.  Our  laser  safety 
guidelines  in  their  present  form  must  be  considered  tentative  and  used 
with  caution  until  the  characteristics  of  the  supersaturation  effect 
are  determined  for  conditions  other  than  those  used  in  our  experiments. 

* 

Present  indications  are  that  c  is  a  universal  constant  for  a 

m 

given  laser  and  path  configuration.  However,  it  might  be  expected  that 
* 

o  would  vary  with  path  height,  since  the  turbulence  structure  changes, 
m 

Also,  laser  variables — such  as  transmitted  beam  diameter,  divergence, 

and  wavelength — might  well  have  an  effect,  as  mentioned  by  Kerr  (1970). 

* 

Mevers  et  al.,  (1969)  have  tentatively  found  that  a  is  less  than  half 

m 

as  large  for  a  wavelength  of  1.15  u  as  it  is  for  0.514  ^ .  This  wave¬ 
length  dependency  should  be  confirmed  and  further  quantified. 


The  observed  variation  of  scintillation  frequency  with  saturation 
should  also  be  examined  in  more  detail,  since  eye  damage  depends  upon 
the  duration  as  well  as  the  intensity  of  a  pulse. 


B 


When  sufficient  data  of  this  nature  are  available,  they  should  be 
incorporated  into  a  more  comprehensive  and  accurate  laser  eye*  safety 
guide.  In  view  of  the  findings  of  our  recent  research,  it  is  probable 
that  previous  guidelines  can  be  safely  revised  to  be  less  restrictive. 


II!  HORIZONTAL  PROPAGATION  EXPERIMENT 


A.  General 

The  specific  objectives  of  this  experiment  were  as  follows: 

(1)  To  verify  the  existence  and  nature  of  the  saturation 
effect  with  regard  to  optical  scintillation,  and  to 
determine  if  the  saturation  level  for  increasing 
thermal  turbulence  is  equivalent  to  that  for  increasing 

i 

range  under  low- turbulence  conditions. 

(2)  To  check  the  quality  of  the  log-normal  fit  to  the 
probability  distribution  of  the  intensity  fluctuations 
under  saturation  conditions. 

(3)  To  establish  firmer  relationships  between  the  meteorolog¬ 
ical  conditions  and  the  strength  of  the  optical 
fluctuations. 

Toward  these  ends,  we  carried  out  an  experimental  program  to 
measure  the  intensity  fluctuations  in  the  horizontal  beam  from  a  3-mtt 
helium/neon  laser  at  a  height  of  1.7  m  above  the  ground.  Two  5-mm- 
aperture  receivers  were  weed,  one  fixed  at  460-n  range  and  the  other 
movable  over  a  range  from  460  to  3500  n.  The  fixod  control  receiver 
sampled  the  beam  simultaneously  with  the  mobile  receiver  by  means  of  a 
small  diagonal  mirror,  and  thus  served  as  a  reference  for  each  set  of 
measurements  at  different  ranges.  In  this  way  the  contr.  receiver 
supplemented  the  independent  meteorological  measurements  in  monitoring 
the  meteorological  conditions  as  the  range  was  varied  for  the  mobile 
receiver  from  run  to  run. 
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B.  Site  Description 


As  shown  in  Figure  2,  the  site  selected  for  this  experiment  is 
located  in  an  area  of  rice  fields  approximately  midway  between  Woodland 
and  Sacramento,  California.  The  site  and  surrounding  area  are  extremely 
level,  with  scattered  low  dikes  and  levees  being  the  only  significant 
relief.  In  addition,  the  ground  cover  is  relatively  uniform  over  much 
of  the  area.  The  particular  field  we  used  is  approximately  2  X  5  kit 
and  had  a  recently  plowed  and  smoothed  bare  surface.  With  this  type  of 
site,  we  could  be  reasonably  confident  that  the  meteorological  conditions 
and  thermal  turbulence  would  be  horizontally  uniform;  hence,  meteorolog¬ 
ical  measurements  at  a  fixed  location  could  be  considered  representative 
of  the  entire  propagation  path.  Some  idea  of  the  flatness  and  uniformity 
of  this  site  may  be  gained  from  the  photograph  in  Figure  3.  The  detailed 
layout  of  the  laser  propagation  path  and  instrumentation  is  illustrated 
in  Figures  4  and  5. 

C.  Instrumentation 

1 .  Laser  Measurements 

A  Spectra-Physics  Model  115  helium/neon  laser  (wavelength 

o 

6328  A),  with  an  output  power  of  approximately  3  mW,  was  used  for 

this  experiment.  The  laser  was  positioned  so  that  the  beam  was  1.7  m 

above  the  ground.  Pointing  was  controlled  by  means  of  a  precision 

table  with  micropositioners .  Suitable  beam-expanding  optics  were  used 

to  obtain  a  transmitted  beam  diameter  of  approximately  30  mm  (between 
_2 

e  “  intensity  points).  Beam  divergence  was  maintained  at  about  0.5  mrad 
(lull  angle),  to  give  a  25-cm  beam  diameter  at  the  control  receiver 
range  of  460  m,  and  a  190-cm  beam  diameter  at  the  farthest  mobile 
receiver  range  of  3500  m.  A  spatial  filter  with  an  aperture  of  20u  was 
used  to  ensure  a  smooth  mean  (Gaussian)  intensity  profile  across  the  beam. 


12 


FIGURE  2  LOCATION  OF  WOODLAND  FIELD  SITE 


FIGURE  3  VIEW  OF  WOODLAND  FIELD  SITE,  LOOKING  NW 


FIGURE  4  PLAN  VIEW  OF  WOOOLAND  FIELD  SITE 
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Two  essentially  identical  receivers  were  used  in  the 
experiment  to  measure  the  signal  simultaneously  at  two  ranges.  As 
shown  in  Figure  6,  the  control  receiver  assembly  was  fitted  with  a 
2-cm-diameter  diagonal  mirror  so  that  the  main  portion  of  the  receiver 
could  be  positioned  outside  of  and  at  right  angles  to  the  beam  propaga¬ 
tion  path.  The  basic  receiver  system  is  illustrated  in  Figure  7.  Each 
receiver  employed  a  5-mm  diameter  aperture  (to  simulate  the  average  size 

o 

of  a  human  eye  pupil),  and  a  10  A  interference  filter  to  reduce  back¬ 
ground  radiation.  An  RCA  Model  7265  photomultiplier  tube  (PMT)  with  an 
S-20  photocathode  surface  was  used  in  each  receiver.  The  output  signal 
from  the  PMT  was  fed  into  a  logarithmic  amplifier,  whose  output  was  then 
recorded  in  EM  mode  on  an  analog  magnetic  tape  recorder.  The  mean  PMT 
output  current  was  monitored  by  a  meter  using  an  R-C  network  with  a  time 
constant  of  approximately  seven  seconds. 

2.  Meteorological  Measurements 

The  primary  meteorological  instrument  was  a  fast-response 
differential  thermometer  (Figure  8)  designed  to  furnish  a  measurement 
of  the  temperature  structure  function.  The  mechanical  design  of  the 
instrument  was  patterned  after  the  differential  thermometer  used  by 
Wright  and  Schutz  (1967).  in  both  instruments  a  vane  is  used  to  main¬ 
tain  the  orientation  of  two  horizontally  separated  temperature  sensors 
perpendicular  to  the  wind  direction.  However,  in  our  device  we  used 
two  0.0005- inch-diameter  chrome 1-constantan  thermocouples,  physically 
separated  by  30  cm  but  electrically  connected  in  series.  The  voltage 
generated  by  the  temperature  difference  between  the  thermocouples 
(64  u  V/°C,  as  determined  from  repeated  calibration  tests)  is  then 
amplified  by  a  high-gain  differential  amplifier  (Tektronix  Model  1A7A) 
and  recorded  in  FM  mode  on  an  analog  magnetic  tape  recorder.  The  time 
constant  of  the  differential  thermometer  was  determined  experimentally 
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FIGURE  6  CONTROL  RECEIVER  (LEFT)  ANO  MOBILE  RECEIVER  SET-UP 
AT  460  m  RANGE 
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FIGURE  7  BLOCK  DIAGRAM  OF  RECEIVER  SYSTEM 


to  be  0.015  s  when  ventilated  at  2  m/s.  This  response  is  fast  enough 
so  that  the  error  in  the  computed  fluctuation  variance  caused  by 
neglecting  higher  frequencies  should  be  small. 

The  differential  thermometer  was  installed  on  a  mast  at  the 
approximate  laser  beam  level  (1.7  m)  near  the  control  receiver,  460  m 
from  the  laser.  Also  located  at  this  point  were  an  anemometer  and  vane 
for  measuring  mean  wind  velocity. 

Extensive  supplementary  micrometeorological  measurements  were 
obtained  under  subcontract  by  the  Department  of  Water  Science  and 
Engineering  of  the  University  of  California  at  Davis.  Their  instrumen¬ 
tation  installation  (Figure  9)  included  Thornthwaithe  anemometers  at 
five  levels  and  aspirated  temperature  sensors  at  nine  levels  to  furnish 
detailed  wind  and  temperature  profiles  up  to  4-m  height.  These  measure¬ 
ments  were  later  used  to  compute  Richardson  numbers  to  serve  as  estimates 
of  atmospheric  stability.  Additional  measurement*  included  soil  tempera¬ 
ture  and  water  vapor  vertical  gradients,  and  net  radiation.  These  were 
used  to  establish  the  surface  heat  budget  and  the  vertical  heat  flux. 
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D. 


Experimental  Procedures 


The  basic  experimental  configuration  is  illustrated  in  Figure  10. 

The  control  receiver  was  located  460  m  from  the  laser,  and  five  positions 
were  established  for  the  mobile  receiver:  460,  1180,  1900,  2700,  and 


FIGURE  10  EXPERIMENTAL  CONFIGURATION  FOR  HORIZONTAL-PATH  EXPERIMENT 

3500  m  from  the  laser.  The  signal  inputs  from  the  two  receivers  and  the 
differential  thermometer  were  brought  into  a  central  recording  van,  which 
was  located  at  about  500  m  range,  near  the  control  receiver  and  the 
meteorological  tower.  The  signal  from  the  mobile  receiver  was  brought 
back  in  FM  mode  from  the  various  receiver  positions  by  means  of  a  3-km 
length  of  field  wire,  with  taps  at  the  appropriate  ranges.  In  this  way, 
dependable  low-noise  data  recordings  were  obtained. 

The  experimental  technique  consisted  of  5-minute  runs  with  both 
receivers  measuring  the  signal  simultaneously  at  different  ranges.  The 
mobile  receiver  was  moved  from  one  point  to  the  next  point  between  runs 
by  means  of  a  light  truck,  which  housed  the  electronics.  After  com¬ 
pleting  a  run,  the  truck  was  driven  to  the  next  point,  the  receiver  was 
aligned  and  connected  to  the  signal  wire,  10  seconds  of  60-Hz  test  signal 
was  transmitted,  and  then  the  recording  of  a  5-min  run,  including  cali¬ 
bration  levels,  was  begun.  This  procedure,  from  the  start  of  one  run  to 
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the  start  of  the  next,  took  approximately  15  to  20  min.  The  runs  where 
the  mobile  receiver  was  at  the  same  range  as  the  control  receiver  served 
to  furnish  a  comparison  of  the  two  outputs,  which  should  be  equivalent 
at  the  same  range. 

Prior  to  each  run,  the  PMT  supply  voltage  on  each  receiver  was 
adjusted  to  give  an  approximate  mean  output  current  of  10  yA.  The  log 
amplifiers  were  adjusted  to  give  zero  output  for  an  input  of  10  |iA,  so 
that  the  means  of  the  recorded  (logged)  signals  would  '->«->  near  zero. 

E.  Data  Analysis 

As  detailed  in  Table  I,  a  total  of  89  five-minute  experimental  runs 
(less  three  aborted  runs)  were  recorded  during  the  period  10-27  November 
1968.  These  data  were  processed  by  both  analog  and  digital  techniques. 
The  primary  aim  was  to  compute  the  standard  deviation  of  the  signal  in 
each  of  the  three  channels — control  receiver,  mobile  receiver,  and 
differential  temperature  (AT) — for  each  of  the  runs. 

A  Meteorology  Research  Incorporated  "Sigma  Meter"  was  used  for  the 
analog  processing.  This  is  a  two-channel  device  that  half-wave  rectifies 
the  incoming  signal,  takes  the  average  of  the  rectified  signal  (which 
represents  |v|,  the  average  signal  deviation  from  the  mean),  and  then 
determines  the  standard  deviation  a  from  the  relatione  =  1.25  |v|, 
which  holds  for  a  Gaussian  distribution.  The  time  constant  for  the 
Sigma  Meter  is  about  30  s  for  the  particular  averaging  time  selected. 

Since  past  work  has  repeatedly  shown  that  measured  log  intensity 
amplitude  probability  distributions  conform  well  to  Gaussian  distribu¬ 
tions,  the  restriction  of  the  Sigma  Meter  on  distribution  shape  did  not 
affect  the  standard  deviations  computed  for  the  two  laser  signal 
channels.  However,  it  was  found  that  the  AT  distributions  were  non- 
Gauss  ian,  and  hence  the  output  oi  the  anal  ig  computer  was  not  a  true 
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Table  I 


USER  DATA  SUMMARY,  WOODLAND  EXPERIMENT 


Date 


11/10/68 

I 

I 

11/11/68 


♦ 

11/12/68 


11/13/68 


11  M/68 


gm 

jpfa 

BUM 

Ree . 

Site 

Weather 

Tape 

No. 

Remarks 

1546 

1 

OlO  160/10 

20 

Foggy  earlier 

1615 

2 

0l0  200/3 

20 

1636 

3 

/0  15 

20 

1331 

1 

100©/©5H  160/10+15 

20 

1350 

1 

Same 

10 

Abort 

1353 

1 

Same 

10 

1412 

2 

100  0511  180/10 

10 

1433 

3 

Same 

10 

1459 

4 

Same 

10 

1524 

5 

100  0511  180/15 

10 

At  small 

1115 

1 

025  315/5 

2 

50 

1133 

2 

Same 

50 

1155 

A 

025  330/10 

50 

1222 

4 

0  25  330/7 

50 

1249 

5 

Same 

50 

1344 

5 

Same 

50 

1404 

4 

0  25  320/7 

50 

mm 

3 

Same 

20 

in 

2 

Same 

20 

1512 

_1_ 

0  2  5  330/10 

20 

Ch.B  low  cal. 

noisy 

1540 

5 

Same 

10 

1327 

1 

Ol5  300/20+25 

3 

20 

1343 

2 

Same 

20 

1403 

o 

w 

/©V®  15  300/20+25 

20 

4 

Same 

10 

5 

Same 

10 

At  small 

1501 

5 

Same 

10 

1532 

4 

/©V®  15  300/15-r20 

10 

At  very  small 

1609 

3 

Same 

10 

1627 

2 

Same 

10 

1644 

1 

/©V  ©  15  290/13+16 

20 

Ch.B  attcn. 

uccid.  moved 

0917 

1 

15O®/015  150/5 

■» 

10 

0930 

2 

150©/©  15  135/8 

10 

09  14 

3 

Same 

10 

UC/Davis  off 

1012 

— 

1500/015  140/6 

10 

1032 

3 

Same 

10 

UC  back  on 

1049 

2 

1500.015  150/7 

10 

1  lr-3 

-i- 

130®  ®  13  150/7 

10 

1122 

2 

i>a  me 

10 

Temp.  48”  ¥ 

i  ■  :s7 

3 

Same 

10 

1152 

4 

.'ime 

10 

1210 

**3~ 

150©  013  1.15/6 

10 

5  min.  extra 

ftn  on  recorder 

1224 

4 

Same 

10 

Li..  B  cals.  10  a 

4  1  O'jf 

1241 

3 

Same 

10 

1254 

•i 

80  ©  7  t  ili  6 

10 

UIOT 

l 

Sa  nv 

10 

Ltshi  ratn 

beginning 
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Table  1  (Continued) 


■ 

Date 

mu 

Brit 

Weather 

(pV/cm) 

AT 

Sens. 

Remarks 

mm 

11/25/68 

1005 

1 

0  20  320/5 

5 

20 

Abort 

1012 

1 

025  300/7 

20 

■9 

1026 

2 

0  20  315/6 

20 

48 

1042 

3 

Same 

20 

Aoort 

49 

1052 

3 

0  20  300/8 

20 

50 

1156 

4 

0  20  300/15 

20 

Run  called 

51  on  tape 

51 

1215 

-5- 

020  300/15+20 

20 

52 

1272 

•1 

020  300/16+22 

20 

53 

1249 

3 

020  310/15+22 

20 

54 

1302 

2 

0  20  310/16 1-22 

20 

55 

1317 

mm 

Same 

20 

Keg.  spikes 

1 

Ch.  B 

56 

1404 

1 

Same 

20 

57 

1429 

_2_ 

020  315/15+20 

10 

AT  small 

58 

11/26/68 

0836 

i 

/Q)20  020/7 

6 

10 

59 

0854 

2 

/Q)20  350/5 

20 

60 

0912 

3 

Same 

20 

61 

0937 

4 

Same 

20 

62 

0955 

-5- 

/-(D)20  045/4 

20 

63 

1011 

4 

/-©20  030/6 

20 

64 

1033 

3 

/-©20  030/5 

20 

65 

2 

/-®20  315/3 

20 

Spikes  Ch.  B 

66 

1158 

-1- 

'-©20  360/10 

20 

67 

1214 

2 

/-®20  330/12 

20 

Spikes  Ch.  B 

68 

1230 

3 

/-©  330/12+16 

20 

69 

1252 

4 

/-©  320/12+16 

20 

70 

4 

/-©  315/10 

20 

Beam  ht.  exp. 

71 

1330 

-5- 

/-©20  320/12 
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?  measure  of  the  standard  deviation.  The  AT  signal  is  quite  intermittent, 

with  quiet  periods  alternating  with  disturbed  periods  as  the  wind  brings 
pockets  of  warm  air  past  the  sensor.  This  gives  rise  to  a  highly  lepto- 
kurtic  frequency  distribution,  as  will  be  illustrated  in  the  next  section. 

The  problem  with  the  AT  distribution  necessitated  reprocessing  the 
data  in  digital  fashion.  Since  it  took  little  additional  effort,  the 
two  laser  signal  channels  were  reprocessed  as  well  as  the  AT  channel. 

The  procedure  involved  digitizing  each  channel  at  a  sampling  rate  of 
approximately  2000  samples/s,  and  then  compiling  the  signal  amplitudes 
into  256-cell  histograms,  from  which  means  and  standard  deviations  could 
be  computed.  In  order  to  reduce  the  influence  of  any  trends  in  the  mean 
levels  upon  the  standard  deviations,  each  run  was  broken  up  into  approxi¬ 
mate  1-min  segments.  The  means  of  each  of  these  segments  were  then  used 
to  correct  for  any  trends,  and  a  composite  distribution  with  a  common 
mean  was  formed  from  the  segment  distributions.  The  standard  deviation 
of  this  composite  distribution  was  then  taken  to  represent  the  entire 
five-minute  run. 

Several  of  the  runs  were  also  processed  using  a  bank  of  third- 
octave  filters  to  obtain  approximate  analog  power  spectral  distributions. 
The  power  spectra  were  used  to  determine  the  dominant  frequencies  in 
the  signals  on  each  of  the  three  channels. 

Each  of  the  signals  was  monitored  in  the  field  on  an  oscilloscope 
(CRO).  It  was  found  desirable  to  take  Polaroid  photographs  of  a  short 
time  sample  of  each  signal  as  displayed  on  the  CRO,  for  later  use  in 
examining  and  editing  the  data.  An  example  of  one  series  of  such  CRO 
time  exposures  is  given  in  Figure  11.  The  width  of  the  signal  envelope 
is  approximately  proportional  to  the  standard  deviation.  It  may  be 
noted  in  the  figure  that  while  the  width  of  the  signal  envelope  for  the 
control  receiver  remains  about  the  same  from  run  to  run,  that  for  the 
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(■)  460  m  —  1346  PST  (mobile  receiver  ediecent  to  (b)  1180  m  —  1322  PST 

control  receiver) 


(c)  1900  m  —  1311  PST  (d)  2700  m  —  1257  PST 


(e)  3600  m  —  1243  PST  ta -7472-20 

FIGURE  11  FIFTEEN-SECOND  OSCILLOSCOPE  TIME  EXPOSURES  OF  LOGGED  SIGNALS 
FROM  MOBILE  RECEIVER  (TOP  TRACE  IN  EACH  PHOTOGRAPH)  AND 
CONTROL  RECEIVER  AT  460  m  (BOTTOM  TRACE),  ON  27  NOVEMBER  1968. 
The  control  receiver  range  is  given  beneath  each  photograph.  Vertical  scale 
14  HB/di vision:  Horizontal  scale  10  ms/division. 


mobile  receiver  decreases  with  increasing  range.  This  feature  of  the 
data  will  be  discussed  in  the  next  section. 

F.  Results 

1.  Comparison  of  Scintillation  Measurements  with  Theory 

The  analysis  procedure  described  in  the  preceding  section 
yielded  values  of  the  standard  deviations  of  the  two  (log-amplif led) 
signals  from  the  control  and  mobile  laser  receivers,  and  of  the  signal 
from  the  differential  thermometer  (AT).  These  values  are  listed  in 
Table  II. 

The  measurements  of  9 (AT)  were  used  to  derive  values  of  C  , 

R 

the  refractive  index  structure  constant,  in  the  following  manner.  The 
temperature  structure  function  is  defined  as 

dt  -  <r  -  t2>2  =  o2(AT)  .  (3) 


The  well  confirmed  Kolmogorov  "two-thirds  law"  predicts  that  in  the 
inertial  subrange. 


D 

T 


2/3 

r 


(4) 


where  C  is  the  temperature  structure  constant  and  r  is  the  distance 
T 

between  measurement  points.  Thus 


C 

T 


a  (AT) 


-1/3 

r 


(5) 


It  can  easily  be  shown  that  the  structure  constants  for 
refractive  index  and  temperature  are  essentially  proportional  and 
related  by 

C  =  (79  x  lo“6)  (p/T2)  C  , 
n  T 


(6) 


28 


Table  II 


RESULTS  FROM  MOQDIAND  DATA  PROCESSING 


Run  No. 

o  (control),  dB* 
m  7 

o  (mobile).  dB 
m 

0<AT).  °c 

1A 

1.19 

1.23 

MM 

IB 

0.80 

1.62 

■1 

1C 

1.20 

3.61 

mSSM 

2 

4.32 

3.45 

USA 

4 

3.48 

3.39 

wBm 

5 

2.61 

4.98 

6 

1.34 

4.26 

■SB 

7 

0.68 

2.33 

8 

0.25 

2.52 

WSSM 

9 

5.30 

5.20 

itsi 

10 

5.37 

5.52 

BUB 

11 

5.20 

4.04 

0.297 

12 

5.29 

3.64 

0.323 

13 

5.13 

** 

0.304 

14 

4.72 

3.44 

0.243 

15 

4.16 

3.44 

0.202 

16 

3.22 

4.59 

0.158 

17 

2.31 

4.87 

0.119 

18 

1.74 

1.97 

0.088 

19 

0.66 

3.83 

0.035 

20 

3.40 

3.67 

0.165 

21 

3.41 

4.81 

0.160 

22 

2.83 

4.63 

0.131 

23 

2.44 

4.77 

0.104 

24 

1.14 

5.14 

0.047 

25 

0.84 

4.14 

0.035 

26 

0.44* 

** 

0.016* 

27 

1.35* 

4.40* 

0.058* 

28 

1.96* 

4.25* 

0.089* 

29 

2.38* 

2.49* 

0.106* 

30 

2.01 

2.24 

0.098 

31 

2.24 

5.28 

0.086 

32 

2.38 

6.06 

0.105 

33 

2.28 

5.18 

0.101 

34 

2.17 

5.49 

0.107 

35 

2.07 

5.21 

0.102 

36 

2.16 

2.39 

o 

c 

37 

1.99 

5.23 

0.093 

38 

1.86 

5.56 

0.077 

39 

1.63 

5.60 

0.074 

40 

1.05 

5.30 

0.053 

41 

o.«o 

2.56 

0.030 

Divide  value*  by  -1.34  te  convert  from  dB  to  conventional 
natural  logarithms 

‘Stable  condition*  (surface-baaed  temperature  inversion) 
“Data  missing  or  invalid 


Tabic  II  (Continued) 


0  (control),  dB 
m 

9  (Mobile),  dB 

o(AT).  °c 

0.51 

1.85 

WBESM 

EH 

0.51 

1.31 

HUB 

0.52 

0.57 

WSSSm 

46 

5.14 

5.46 

mSSM 

47 

4.73 

5.16 

BS3H 

49 

5.38 

4.60 

SO 

4.83 

4.30 

— 

51 

4.61 

3.55 

B2H 

S2 

4.37 

4.00 

HBH 

S3 

4.15 

4.49 

hbh 

54 

4.12 

4.85 

W 

55 

4.03 

2.89 

'SIS 

56 

2.74 

2.85 

0.112 

57 

1.91 

** 

58 

1.62 

1.75 

59 

2.49 

5.33 

sii 

60 

2.85 

5.74 

61 

3.18 

5.07 

62 

4.41 

3.50 

63 

4.44 

4.06 

0.204 

64 

3.60 

5.30 

0.190 

65 

4.29 

** 

0.211 

66 

4.40 

4.50 

0.211 

67 

4.27 

** 

0.206 

68 

4.23 

4.95 

0.203 

69 

2.74 

5.38 

0.122 

70 

2.02 

5.27 

0.085 

71 

2.24 

5.13 

0.102 

72 

1.23 

5.02 

0.051 

73 

1.29 

4.42 

0.044 

74 

0.79 

2.01 

0.028 

75 

0.86 

0.92 

0.029 

76 

4.96 

5.37 

0.217 

77 

5.46 

5.69 

0.199 

78 

5.37 

4.78 

0.268 

79 

5.48 

3.70 

0.251 

80 

3. 18 

2.80 

0.288 

81 

3.40 

4.08 

0.274 

82 

3.23 

4.61 

0.268 

83 

4.58 

5.06 

0.237 

84 

4.84 

5.00 

0.218 

83 

4.39 

3.82 

0.217 

86 

4.59 

4.97 

0. 190 

87 

3.93 

4.81 

0. 183 

HM 

3.39 

4.53 

0.135 

m* 

0.86 

•  • 

0.046 
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where  p  is'-atmospheric  pressure  in  millibars  and  T  is  absolute  air 
temperature  in^^kv  Combining  Eqs.  (5)  and  (6),  we  have 


C  =  (79  X  l(f6)  (p/T2)  (r"1/3>  o(AT) 
n 


(7)  / 


Based  upon  a  value  for  r  of  30  cm  for  our  differential  thermometer’ 

✓ 

which  is  the  distance  between  the  thermocouple  probes,  and  atypical 

y-' 

value  of  1000  millibars  for  p,  we  obtain 


C  =  0.116  a(AT)  /T  ,  (8) 

n 

a  o  -1/3 

where  a (At)  and  T  are  in  units  of  K  and  C  is  in  units  of  m  .  This 

n 

is  the  equation  that  was  used  to  calculate  C  . 

n 

The  derived  values  of  C  are  plotted  versus  the  measured  log 

n 

intensity  standard  deviations  at  the  control  receiver,  a  (control),  in 

m 

Figure  12.  The  control  receiver  was  always  at  a  fixed  range  of  460  m. 

The  dashed  line  in  the  figure  is  the  least-squares  fit  to  the  data  for 
-6  -1/3 

C  <  0.35  X  10  m  .  Values  for  higher  C  show  evidence  of  saturation, 
n  n 

even  at  the  460  m  range,  and  hence  were  excluded  from  the  fit,  although 
they  are  shown  on  the  plot.  Knowing  that  the  range  R  is  460  m,  and  that 
the  wavelength  X  is  0.6328u,  the  slope  of  the  fitted  line  can  be  used  to 
derive  the  coefficient  in  Tatarski's  theoretical  expression  for  log-normal 
scintillation,  o^  (Eq.  1).  The  fitted  form  of  the  equation  comes  out  to 


a  .  1.06  „U/ia  («/,., 7/13  C 
t  n 


Tatarski  gives  theoretical  values  that  yield  comparable  coefficients  of 
1.11  for  a  plane  wave  and  0.72  for  a  spherical  wave.  The  measured  value 
of  1.06  agrees  well  with  these  values. 
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TA-7472-26 

FIGURE  12  MEASURED  RELATION  BETWEEN  THE  LOG  INTENSITY 
STANDARD  DEVIATION  AT  THE  CONTROL  RECEIVER 
AT  460  m  RANGE,  AND  THE  REFRACTIVE  INDEX 
STRUCTURE  CONSTANT  (Cn)  DERIVED  FROM  THE 
DIFFERENTIAL  THERMOMETER  DATA 


The  striking  feature  of  Figure  12  is  the  good  fit  obtained  to 

Tutarski's  equation  for  short  range  and  sraal 1-to-moderate  values  of  C  . 

n 

The  scatter  of  the  data  points  is  fairly  small,  considering  that  (1)  the 

measurements  for  a  (control)  and  C  were  obtained  by  completely  indepen- 
m  n 

dent  optical  and  meteorological  techniques,  and  (2)  the  laser  measurements 
represent  a  line  integral,  while  the  C^  values  are  based  upon  meteorolo¬ 
gical  measuremenss  at  a  single  point, 

l' sing  Kq.  (9),  values  of  were  calculated  for  each  value  of 
measured  by  the  mobile  receiver  at  the  various  ranges  (160,  1180, 

Kt 

190o,  2700,  and  OatH)  ra),  and  plotted  versus  0  (mobile)  in  Figure  13. 

IQ 

The  data  cUarlv  show  not  only  that  there  is  a  saturation  of  scintillation 
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<T  (THEORETICAL) 


TA-7472-27 


FIGURE  13  COMPARISON  OF  MEASURED  VALUES  OF  LOG 

INTENSITY  STANDARD  DEVIATION  AT  THE  MOBILE 
RECEIVER  WITH  VALUES  PREDICTED  BY  TATARSKI'S 
THEORY,  SHOWING  THE  SUPERSATURATION  EFFECT 


at  about  =  2,  but  also  that  there  is  a  supersaturation  region 

beyond,  in  which  a  decreases  with  increasing  a  .  The  dashed  curve  in 
m  —————  t 

the  figure  represents  the  best  least-squares  fit  obtained  with  the  model 

a 


a  - 
m 


1  +  aa 


(10) 


The  equation  for  the  fitted  curve  comes  out  to  be 


2 

1  +  0.16  rj 


(ID 
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* 

This  function  has  a  maximum  scintillation  magnitude  of  o  =  1.2?  at 

in 

O  =  2.55.  In  the  limit  fo*-  small  0,0  =  (0  )  and  the  scirtilla- 

t  t'  n  t  small 

tion  is  represented  by  the  Tatarski  expression.  In  the  limit  for  large 


o  ,  we  have 
t' 


o 

m 


=  6. 49/(0  V 

t  large 


(12) 


Figure  14  shows  the  result  when  the  curve  fitted  to  our  measure¬ 
ments  (Eq.  11)  is  compared  with  the  observations  by  Gracheva  (1967).  The 
function  is  higher  than  the  Russian  measurements  in  the  vicinity  of  the 
peak  scintillation  values,  and  slightly  lower  for  large  The  reason 


FIGURE  14  MEASURED  (om)  VERSUS  THEORETICAL  (o,)  VALUES  OF  LOG 
INTENSITY  STANOARO  DEVIATION.  AS  MEASURED  BY 
GRACHEVA  (1967)  AND  CORRECTED  BY  Oc  WOLF  (1968).  The 
curvet  represent  wtcut  theoretical  tormuletion*  by  (A)  Teienki  (1961). 
(8)  Tetanki  (1966).  end  <C>  De  WoH  (1968).  pfcrt  the  empirical 
(unction  luted  to  the  SRI  date  ID). 
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for  this  disparity  is  unknown,  but  it  is  possible  that  the  Russian 

equipment  lacked  sufficient  dynamic  range  to  measure  large  o  ,  since 

m 

they  did  not  use  log  amplifiers.  The  Russian  experiment  also  employed 
a  mercury  vapor  lamp  rather  than  a  laser. 

In  addition  to  the  Russian  results  (which,  incidentally,  the 
experimenter  did  not  recognize  as  showing  mpersaturation),  other 
recent  experimental  work  by  Deitz  and  Wright  (1969),  Movers  et  al.  (1969), 
and  Ochs  (1969)  also  show  indications  of  the  suporsaturation  effect. 

The  supersaturation  effect  perhaps  can  be  visualized  more 
clearly  by  referring  to  Figure  15,  in  which  the  measured  values  of  a 

m 

are  grouped  according  to  various  class  intervals  of  C  and  plotted  as 

n 

a  function  of  range.  The  scatter  of  the  data  points  in  these  groups  is 

primarily  due  to  the  widths  of  the  class  intervals  used  for  C  .  It  may 

n 

♦ 

be  noted  that  the  value  of  maximum  scintillation  a  is  approximately 

m 

equal  to  1.27  and  does  not  depend  upon  C  .  However,  the  range  at  which 

n 

the  peak  scintillation  occurs  is  a  function  of  C  ;  for  larger  C  the 

n  n 

range  to  the  peak  decreases.  A  surprising  result  is  that  at  a  given 

range,  say  3  km,  the  scintillation  will  generally  be  greater  for  small 

C  than  for  large  C  ;  it  will  reach  the  maximum  value  at  some  intermediate 
n  n 

value  of  C  . 

n 

It  should  be  noted  that  the  dashed  curves  in  the  fou;  graphs 
in  Figure  15  represent  the  same  analytical  function  derived  from  all 
the  data,  rather  than  individual  fits  to  the  classed  data. 

The  results  shown  in  Figures  13  and  15  are  quite  significant 

in  the  eye-safety  context.  The  maximum  eve-damage  probability  is 

* 

obtained  at  the  peak  value  of  o  -  1.27,  and  at  long  ranges  the  eye- 

18 

damage  probability  goes  down,  rather  than  increasing  approximately 
linearly  with  range,  as  Tatarski's  equation  predicts.  The  new  results 
are  incorporated  into  the  procedure  for  estimating  eye-damage  probability, 
which  is  discussed  in  Section  V. 


35 


FIGURE  15  MEASURED  LOG  INTENSITY  STANDARD  DEVIATION  AT  THE  MOBILE  RECEIVER  AS  A  FUNCTION  OF  RANGE,  FOR 
VARIOUS  CLASS  INTERVALS  OF  Cp.  The  dashed  curve  in  each  graph  represents  the  empirical  function  a  =  0/CI+O.I60 


2.  Log-Normality  ol'  Probability  Distributions 


For  the  eye-damage  estimates  we  assume  that  the  intensity 
fluctuations  follow  a  log-normal  probability  distribution,  based  upon 
the  results  of  our  first,  year's  work.  In  general,  the  validity  of  tins 
assumption  is  confirmed  by  the  new  results.  To  illustrate  this,  a 
reasonably  representative  run  has  been  selected  and  is  depicted  in 
Figures  16-18.  The  probability  distributions  for  the  individual  one- 
minute  segments  of  the  five-minute  run  are  graphed  in  Figure  16  for  the 
control  and  mobile  receivers.  The  composite  five-minute  probability 
distribu lions  for  the  two  receivers  are  compared  with  Gaussian  distribu¬ 
tions  with  the  same  standard  deviations  in  Figure  17,  and  are  compared 
with  each  other  in  Figure  18. 

As  indicated  in  Figure  17,  the  actual  probability  distributions 
of  the  log  intensity  for  the  two  receivers  are  fair  approximations  of 
the  Gaussian  curve.  This  is  the  case  even  though  supersaturation  is 
clearly  occurring  in  the  signal  at  the  mobile  receiver  at  3500  m,  as 
shown  in  Figure  18  by  the  smaller  standard  deviation  for  this  distribu¬ 
tion.  It  may  be  noted  in  Figure  17  that  the  distribution  for  the  mobile 
receiver  has  a  slight  negative  skew.  This  was  found  to  be  typical  for 
detectors  that  were  well  centered  in  the  laser  beam.  (Our  laser  had  a 
Gaussian-shaped  beam  intensity  profile.)  When  the  detector  was  off- 
center  ’ey  a  moderate  amount  with  respect  to  the  beam,  the  probability 
distribution  becomes  positively  skewed,  as  illustrated  by  the  example  in 
Figure  19.  Here  the  bean,  was  apparently  slowly  drifting,  such  that  the 
originally  well  centered  detector  ended  up  near  the  edge  of  the  beam  at 
the  end  of  the  five-minute  run.  The  one-minute  probability  distributions 
start  out  fairly  symmetrical  and  Gaussian,  but  become  increasingly  posi¬ 
tively  skewed,  probably  duo  to  the  proximity  of  the  lower  bound  of  zero 
intensity.  Apparently  a  true  log-normal  distribution  occurs  only  when 
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FIGURE  16  ONE-MINUTE  PROBABILITY  DISTRIBUTIONS  FOR  (a)  CONTROL 
RECEIVER,  ANO  <b)  MOBILE  RECEIVER  AT  3500  m  RANGE, 
FOR  RUN  51. 
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FIGURE  17  COMPOSITE  FIVE-MINUTE  °ROBABILITY  DISTRIBUTION 
AND  GAUSSIAN  FIT  FOR  (a)  CONTROL  RECEIVER,  AND 
(b)  MOBILE  RECEIVER  AT  3500  m  RANGE,  FOR  RUN  51 
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LOG  RELATIVE  INTENSITY  -  06 

TA-7472-31 

FIGURE  18  COMPARISON  OF  COMPOSITE  PROBABILITY  DISTRIBUTIONS 
FOR  CONTROL  AND  MOBILE  RECEIVERS  FOR  RUN  51. 


the  detector  is  located  slightly  off-center  with  respect  to  the  beam 

axis. 

The  standard  deviation  of  the  distributions  shown  in  Figure  19 
might  be  expected  to  increase  as  the  detector  moves  off-center.  As 
shown  in  Table  III,  there  is  some  indication  that  this  may  occur,  but 
the  increase  in  scintillation  is  apparently  sufficiently  small  that  this 
effect  can  be  safely  ignored  in  our  analysis. 


Figure  20  illustrates  the  characteristic  non-Gaussian,  lepto- 
kurt  ie  nature  of  the  probability  distributions  for  At.  This  is  caused 
by  the  fact  that  t  lie  AT  signal  is  typically  quite  intermittent,  with 
periods  of  very  small  fluctuations  interspersed  with  intervals  of  rela¬ 
tively  large  fluctuations.  It  turns  out  that,  as  shown  in  the  example 
in  Figure  21,  the  observed  characteristic  A'f  distribution  is  fitted 
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FIGURE  19  EXAMPLE  OF  CHANGE  IN  SHAPE  OF  PROBABILITY 
DISTRIBUTIONS  AS  BEAM  DRIFT  CAUSES  DETECTOR 
LOCATION  TO  CHANGE  TO  EDGE  OF  BEAM  (RUN  12. 
MOBILE  RECEIVER,  RANGE  2700  m). 


Table  III 


MEANS  AND  STANDARD  DEVIATIONS 
FOR  RUN  12  (MOBILE  RECEIVER) 


Interval 

No. 

Mean  Signal 
Level  (dB) 

Standard 

Deviation  (dB) 

1 

-2.07 

4.36 

2 

-2.37 

4.28 

3 

-3.44 

4.66 

4 

-3.62 

4.45 

5 

-4 . 68 

4.57 

T1 


RELATIVE  DIFFERENTIAL  TEMPERATURE  —  DEG  C 
(•I 


TA-747I-33 

FIGURE  20  (a)  ONE-MINUTE  AND  (b)  COMPOSITE  FIVE-MINUTE 

PROBABILITY  DISTRIBUTIONS  FOR  AT  FOR  RUN  51 
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RELATIVE  PER  CENT  OCCURRENCE  RELATIVE  PER  CENT  OCCURRENCE 


-1.5  -1.0  -0.5  0  0.5  1.0 


DIFFERENTIAL  TEMPERATURE  —  DEC  C 


FIGURl  21  EXAMPLE  OF  FIT  OF  COMPOUND  NORMAL  DISTRIBUTION 
TO  OBSERVED  IT  PROBABILITY  DISTRIBUTION.  (»)  Fitted 
compound  normal  distribution;  (b)  Components  of  fitted  compound 
normal  distributions.  The  parameters  n,  and  n?  represent  the 
total  number  of  occurrences  in  each  of  the  component  distributions 
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reasonably  well  by  the  sum  of  two  normal  distributions  with  different 
standard  deviations. 

3.  Spectra  of  Scintillation 

The  partition  of  laser  signal  fluctuation  energy  in  terms  of 
frequency  was  examined  for  several  sets  of  runs  by  processing  the  analog 
data  through  a  bank  of  third-octave  filters.  This  procedure  yielded 
values  representing  relative  log  intensity  standard  deviation  for  eight 
frequency  bands  from  8  to  1000  Hz.  These  approximate  power  spectra  are 
presented  in  Figure  22  for  two  sets  of  runs  obtained  during  reasonably 
stationary  meteorological  conditions.  Figure  22(a)  represents  a  light- 
wind  case,  when  the  cross-path  wind  component  was  only  about  1  m/s.  The 
spectra  for  the  signal  from  the  control  receiver  at  460  m  are  very  con¬ 
sistent  from  run  to  run,  for  the  first  three  runs  of  the  series,  with 
the  peak  energy  at  about  16  Hz.  The  signal  spectrum  for  the  fourth  run 
(#83),  although  similar,  probably  shows  the  effect  of  an  increase  in 
cross -path  wind  component.  However,  the  mobile  receiver  signal  spectra 
show  the  effects  of  supersaturation  as  the  receiver  is  moved  downrange. 

With  increasing  range,  the  total  fluctuation  energy  decreases,  and  the 
dominant  fluctuations  shift  to  higher  frequencies,  which  was  apparent 
from  inspection  of  the  signals  while  the  data  were  being  obtained.  The 
spectra  for  the  mobile  receiver  also  become  bimodal,  with  two  peaks,  one 
at  about  16  Hz  and  the  other  at  about  150-200  Hz. 

Figure  22(b)  shows  the  power  spectra  for  a  higher-wind  situation, 
in  which  the  cross-path  component  was  about  4  m/s.  The  control  receiver 
spectra  all  peak  at  about  123  liz,  but  the  spectra  for  the  mobile  receiver 
are  quit,  similar  to  those  for  the  light-wind  case  [Figure  22(a)}.  This 
indicates  that  the  wind  speed  has  little  effect  upon  fluctuation  fre¬ 
quencies  in  the  Mipersatur.it  ton  region,  although  it  is  significant  in 
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APPROXIMATE  ANALOG  POWER  SPECTRA  FOR  CROSS-PATH  WIND 
COMPONENT  OF  I#!  I  *»-9«c.  4NO  M  I  «K 


the  subsaturation  region,  as  Tatarski  (1961)  predicted  and  Ryznar  (1965) 
experimentally  confirmed.  The  reason  for  the  double  peak  in  the  spectra 
for  the  signals  in  the  supersaturation  region  is  unknown,  but  the  lower 
frequency  peak  could  be  associated  with  beam  wander.  It  is  possible 
that  these  spectra  may  oifer  a  clue  into  the  mechanism  for  the  super¬ 
saturation  effect.  This  phenomenon  is  probably  associated  with  multiple 
turbulent  scattering,  in  which  the  high-and-low- intensity  regions  in  the 
beam  are  eroded  by  the  process  of  destructive  interference.  Dcitz  and 
Wright  (1969)  have  offered  some  qualitative  suggestions  along  this  line. 

•1.  Comparison  of  Standard  Deviations  Obtained 
by  Digital  ar.d  Analog  Processing 

The  results  of  the  digital  data  processing  procedure  for 
determining  standard  deviations  were  used  exclusively  in  the  Woodland 
data  analysis  for  deriving  the  results  just  presented.  However,  it  is 
worthwhile  determining  the  degree  of  correlation  between  the  standard 
deviations  derived  by  the  digital  and  analog  methods.  Comparisons  of 
J  (digital)  and  z  (analog)  are  presented  in  Figure  23  The  indicated 
lunar  least-squares  fits  give  the  following  relationships: 

|:(control)J  0.087  -  1.01  (a (control)]  (13) 

dig  anal 

l :  (mobile))  0.082  -  1.02  Immobile)]  (14) 

dig  ana'. 

{'t:.T)i  0.06  •  1.37  (7(aT)l  .  (13) 

dig  anal 

S  m e  tin-  above  relationships  lit  the  data  well,  with  rela- 
!v  ii?tle  i.tt.r,  thi-y  haw  proven  use  fill  tr  alibi .ling  the  .Mtgna 
V. »  r  to  lorreet  the  analog  data  t i or  the  control  rectiver  and  dit- 
u!  t hertnor.etor  collected  d./.ng  the  slant-path  experiment. 
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IV  EXPERIMENT  B— SLANT-PATH  PROPAGATION  EXPERIMENT 


A.  General 

The  specific  objectives  of  this  experiment  were  as  follows: 

(1)  To  verify  the  existence  of  scintillation  magnitude 
saturation  for  slant-path  laser  propagation,  and  to 

* 

examine  the  magnitude  of  the  scintillation  maximum  o 

m 

and  its  dependence  upon  path  geometry. 

(2)  To  develop  relationships  by  means  of  concurrent  measure¬ 
ments  among  slant-path  scintillation,  scintillation  on  a 
horizontal  path  near  the  ground,  and  atmospheric 
conditions. 

(3)  To  determine  the  effect,  if  any,  of  the  pseudo-wind  due 
to  aircraft  motion  upon  scintillation  magnitude. 

To  accomplish  these  aims,  we  fielded  an  experimental  program  at 
Brooks  Air  Force  Base,  Texas,  during  the  period  14-23  April  1968.  The 
main  portion  of  the  experiment  involved  tracking  an  airborne  receiver 
mounted  on  an  orbiting  C-97  airplane  with  a  ground-based  laser. 

B.  Site  Description 

A  plan  view  of  the  area  surrounding  the  experimental  site  is  pre¬ 
sented  xn  Figure  24.  The  average  elevation  is  about  575  ft  above  sea 
level.  The  terrain  is  generally  flat  to  slightly  rolling  in  the  area 
south  and  east  of  the  ground  site  where  the  flights  were  conducted.  The 
vegetation  consists  of  patches  of  scrub  brush  and  low  trees  intermixed 
with  frequent  open  grassy  fields  and  a  few  surbuban  housing  developments. 
The  central  business  district  of  the  city  of  San  Antonio  lies  7.5  miles 
to  the  northwest  of  the  ground  site. 
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FIGURE  24  EXPERIMENTAL  SITE  LOCATION  ON  BROOKS  AFB.  Three  aircraft  ground  tracks  used  in  the 
experiment  are  shown. 


The  immediate  area  around  the  ground  site  was  level  and  covered 
with  short  grass.  As  shown  in  Figure  25,  the  recording  van  was  situated 


FIGURE  23  SITE  LAYOUT  FOR  SLANT-PATH  EXPERIMENT  AT  BROOKS  AIR  FORCE  BASE 


near  the  laser  tracking  mount,  as  was  the  meteorological  instrumentation 
ar.d  the  receiver  for  the  ground  laser. 

C  Instrumentation 

1.  laser  Measurements 

Two  Spectra-Phy  *'.cs  Model  122  3-mW  helium/neon  lasers  wvre 
used  in  this  experiment,  one  for  the  slant-path  propagation  link  and  the 
other  for  a  350-m  horizontal  control  link.  The  instrumentation  setup 
for  the  horizontal  control  path  at  about  1.7  ■  above  the  ground  was 
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essentially  identical  to  that  used  in  the  Woodland  experiment,  again  with 
a  25-cm-diameter  laser  beam  impinging  upon  a  receiver  with  a  5-mm  aperture. 
However,  the  slant-path  portion  of  the  experiment  required  the  development 
of  special  instrumentation. 

A  surplus  antiaircraft  gun  director  mount  was  precisionized  and 
fitted  with  an  azimuth  motor  drive  and  optical  sighting  devices  to  permit 
its  use  in  tracking  the  receiver  on  the  orbiting  aircraft  with  the  laser 
(Figure  26).  The  same  model  laser  was  use.>d  here  as  was  used  in  the  hori¬ 
zontal  control  link.  The  tracking  mount  proved  to  be  capable  of  consis¬ 
tently  accurate  tracking  to  within  about  plus  or  minus  one-half  plane- 
length  (approximately  ±  5  mrad)  with  an  experienced  operator  when  the 
airplane  was  in  smooth,  level,  orbiting  flight.  In  consideration  of  this 
tracking  accuracy,  the  laser  beam  was  diverged  to  approximately  10  mrad. 
The  maximum  azimuthal  tracking  rate  of  the  mount  was  approximately  6°/s. 

The  left  front  side  window  of  a  USAF  C-97  aircraft  was  removed 
and  the  laser  receiver  assembly  installed  in  the  opening,  as  shown  in 
Figure  27.  In  order  to  reduce  the  background  noise  level  in  the  re- 
ceivor-Teausod^by  reflected  sunlight,  it  was  necessary  to  minimize  the 
receiver  field  of  view.  A  value  of  10  mrad  was  found  to  be  suitable. 

This  required  that  the  receiver  assembly  (Figure  28)  be  mounted  on  a 
pivot  with  appropriate  sighting  devices,  so  that  it  could  be  manually 
aimed  at  the  laser  on  the  ground.  A  50- mm  comer-cube  reflector  was 
mounted  on  the  front  portion  of  the  receiver  assembly  to  provide  a 
visual  indication  of  tracking  accuracy  *o  the  ground  tracking  mount 
operator.  An  audio  monitor  of  the  received  signal  level  was  used  by  the 
airborne  receiver  operator  to  aid  in  alsing  the  receiver.  Figure  29 
shows  the  airborne  receiver  and  associated  electronics  installation  in 
the  aircraft,  as  well  as  the  receiver  assembly  in  its  operating  position. 
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FIGURE  28  LASER  TRACKING  MOUNT 
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FIGURE  27  to)  BROAD  VIEW.  AND  ft>)  CLOSE-UP 
VIEW  Of  LASER  SIGNAL  RECEIVER 
INSTALLATION  IN  W  NDOW  OPENING 
OF  U.S.  AIR  FORCE  C-9?  AIRCRAFT. 
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FIGURE  28  AIRBORNE  RECEIVER  AIMING  ASSEMBLY 


Provision  was  made  for  automatic  sequencing  of  a  test  signal 
and  the  calibration  levels  at  the  beginning  and  end  of  each  run,  which 
was  controlled  by  buttons  in  the  hand  grip  or.  the  receiver.  The  basic 
receiver  electronics  system  was  the  same  as  that  used  in  the  Wcodland 
experiment,  with  the  signs',  from  an  RCA  7265  PUT  (with  5-mn-aperture 
diane.  r)  fed  into  a  logarithmic  amplifier.  The  output  of  the  log 
amplifier  was  converted  to  FM  mode  and  transmitted  back  to  the  recording 
van  at  the  ground  by  ueans  of  a  telemetry  link,  which  proved  to  be  very 
trouble-free.  Communications  between  the  ground  station  and  the  airplane 
were  carried  out  using  a  UHF  channel. 
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FIGURE  29  U)  AIRBORNE  RECEIVER  AND  ELEC¬ 
TRONICS  INSTALLATION  INSIDE  C-9? 
AIRCRAFT.  AN*>  tbl  RECEIVER  AIMING 
ASSEMBLY  •**  OPERATING  POSITION 


2.  Meteorological  Measurements 


Ground  instrumentation  used  in  this  experiment  for  monitoring 
the  meteorological  conditions  consisted  of  an  anemometer  and  wind  vane 
for  measurinf  mean  wind  velocity,  and  the  differencial  thermometer  pre- 
viously  described  for  measuring  the  temperature  structure  function.  All 
instruments  were  positioned  at  the  mean  laser  beam  height,  1.7  m. 

Mean  temperature  profiles  up  to  8000  ft  MSL  were  obtained  by 
means  of  a  thermistor  thermometer  installed  on  the  097  aircraft.  The 
shielded  sensor  tube  was  exposed  to  the  airstream  by  a  mount  on  a  35-cm- 
long  boom  attached  to  the  window  clamp  for  the  receiver  assembly  (Figure 
28).  Temperature  profiles  were  recorded  during  spiral  ascents  and  de¬ 
scents  at  the  start  and  end  of  each  flight  by  means  of  a  chart  recorder 
[shown  at  right  bottom  of  Figure  29(a)]. 

D.  Experimental  Procedures 

The  basic  configuration  for  the  slant-path  experiment  is  illustrated 
in  Figure  30.  The  experimental  procedure  consisted  of  conducting  a  series 
of  runs  during  which  the  airborne  receiver  was  tracked  with  the  laser  as 
the  C-97  flew  an  almost  continuous  circular  path  in  a  counterclockwise 
direction  around  the  ground  station.  The  tracking  occurred  over  an 
approximate  180°  arc  from  southwest  to  northeast  of  the  ground  station. 
Several  aircraft  altitudes  between  1500  and  8000  ft  above  ground  were 
employed,  at  various  ground  distances  between  0.5  and  2.5  miles  from  the 
laser.  The  runs  varied  between  2  and  7  minutes  in  duration. 

The  telemetered  signal  from  each  run  was  recorded  in  FM  mode  on  a 
30  in/s  analog  magnetic  tape  recorder  with  a  10  kHx  frequency  bandwidth. 

The  signal  from  the  receiver  cm  the  horizontal  control  path  and 
the  &T  signal  were  recorded  continuously  in  FM  mode  on  a  separate  7-12 
in/s  magnetic  tape  recorder.  In  addition,  the  signals  were  also  ted  into 
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FIGURE  30  EXPERIMENTAL  CONFIGURATION  FOR  SLANT-PATH  MEASUREMENTS 


the  two  channels  of  the  Sigma  Meter,  whose  output  was  displayed  on  a  chart 
recorder.  This  provided  a  real-time  indicator  of  turbulence  conditions. 

A  chronological  summary  of  activities  during  the  slant-path  experi¬ 
ment  is  presented  in  Table  IV.  Five  days  were  required  for  the  equipment 
installation  and  testing,  followed  by  two  days  of  data  collection.  The 
quantity  of  data  collected  during  the  field  program  was  limited  by  the 
restricted  availability  of  the  C-97  aircraft. 

E.  Data  Analysis 

A  fairly  substantial  number  of  runs  were  completed  on  the  two  days 
of  19-20  April  1969,  as  summarized  in  Table  V.  The  slant-path  signal 
data  for  each  run  wore  digitized  at  a  rato  of  approximately  15,000 
samples/s  and  converted  into  amplitude  probability  distributions 
(histograms)  essentially  as  was  done  with  the  Woodland  data.  However, 
only  one  histogram  per  run  was  computed  for  the  slant-path  data. 


* 
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Table  IV 


CHRONOLOGICAL  SUMMARY,  SLANT- PATH  EXPERIMENT 


Day 

Activity 

14  April 

Unpacked  equipment;  began  ground  installation 

15  April 

Completed  ground  installation 

16  April 

Completed  C-97  aircraft  installation 

17  April 

Test  flight,  operator  tracking  practice 

18  April 

Test  flight,  6  preliminary  data  runs  collected 

19  April 

Data  flight,  4  hours,  24  runs  recorded 

20  April 

Data  flight,  2-1/2  hours.  ?  uns  recorded 

21  April 

Aircraft  in  maintenance 

22  April 

Aircraft  in  maintenance 

23  April 

Data  taken  on  ground  link;  aircraft  grounded 
because  of  mechanical  problems  and  unavailable 
remainder  of  week  because  of  other  commitments 

The  principal  problem  in  the  analysis  of  these  data  resulted  from 
the  intermittent  nature  of  the  signal  as  the  tracking  was  on  and  off 
target.  When  the  tracking  was  off  target,  the  signal  consisted  only  of 
background  noise.  This  can  be  seen  in  Figure  31,  which  shows  examples 
of  15-s  time-exposure  photographs  of  the  CRO  display  for  selected  runs. 
At  the  longer  ranges,  the  signal-to-noise  ratio  (SNR)  was  frequently  low 
and  the  laser  signal  onvelope  merged  with  the  background  noise  envelope. 
These  two  distributions  are  superimposed  upon  one  another,  as  indicated 
in  the  sample  histograms  shown  in  Figure  32.  In  the  first  example  in 
this  figure  the  signal  and  noise  distributions  are  rather  clearly  sepa¬ 
rated,  but  there  were  many  runs  similar  to  the  second  example  where  the 
w-  distributions  were  merged  together  more  closely. 
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Table  V 


1ASKK  DATA  SUMMARY,  SLANT-PATH  EXPERIMENT 


HHi 

Da  l  !• 

Kmul 

Al  re  raft 

Altitude 

(ft  above 
round ) 

Ground 

Distance 

(ml  1cm) 

giffl 

SslB 

P 

Nuabera 

Weather,  Remarks 

Bfflni 
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n 

1/ 19/  G9 

1018: 00 

5000 

n 
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4 

3 

O  35  130V9  nph  78"  F 

El 

1026;  15 

3000 

is 

3000 

4 

3 

Same 

A 

1036:15 
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mSM 
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4 

3 
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•1 

10-17*  13 
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MEM 

*  -ill 

4 

3 

Same 

5 

1057;  15 
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KB 

5 

3 
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G 

1107: 10 
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In 

l  fell 

5 

3 

O  25  180/10  70° 

7 

1120:  15 ' 

3000 

mem 

1850 

3 

3 

S-IB0 

H 

i  ’.29:  00 

3000 

1.0 

1650 

5 

3 

So  aw 

9 

1 136: 40 

J000 

1.0 

1850 

5 

3 

O  »3  135/12 

10 

11-17:00 

1500 

950 

5 

4 

Ground  opor.  couldn't  track 

11 

1155: 15 

3000 

1.0 

1950 

6 

4 

O  I*  140/10 

12 

1206:  IS 

6U00 

MSM 

3150 

6 

4 

Same 

U 

1331*00 

GOoO 

2450 

6 

1 

O  25  130/13.17 

M 

1342:  30 
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u 

3150 

6 

4 

Same 

13 

1352:00 

4000 

19 

2400 

mm 

4 
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16 

1106:05 
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2400 

| 

4 

Abort 

17 

1410:00 
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19 
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4 
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18 
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m 
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B 

4 
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iy 
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mM 

4 
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20 
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9 

5 
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Mm 
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22 
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SI 
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8 

5 
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1350 

8 

5 
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9 
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FIGURE  32  SAMPLE  HISTOGRAMS  FOR  SLANT-PATH  EXPERIMENT 
SHOWING  (a)  GOOD  SIGNAL-TO-NOISE  RATIO  (SNR), 
AND  (b)  POOR  SNR 
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The  analysis  technique  used  to  find  the  standard  deviation  of  the 
signal  distribution  was  simply  to  locate  the  mode  (peak  value)  of  this 
distribution  on  the  histogram— if  there  was  a  distinguishable  mode— and 
then  to  use  the  high-amplitude  half  of  the  distribution  to  graphically 
compute  the  standard  deviation,  assuming  that  the  distribution  is  Gaussian. 
It  is  recognized  that  this  procedure  will  give  only  approximate  results. 
More  accurate  values  should  be  obtained  when  a  new  technique  currently 
under  development  is  completed.  This  new  procedure  involves  the  use  of 
a  nonlinear  least-squares  computer  program  to  fit  a  model  distribution 
to  the  entire  observed  histogram.  One  of  the  model  parameters  obtained 
from  this  fitting  routine  is  the  standard  deviation  of  the  signal  distri¬ 
bution.  The  model  attempts  to  take  into  account  and  remove  the  influence 
of  the  "artificial"  signal  intensity  fluctuation  that  occurs  when  the 
beam — which  has  a  Gaussian  intensity  profile — is  scanned  across  the  re¬ 
ceiver  in  the  course  of  normal  tracking. 

The  chart  records  of  the  standard  deviations  produced  by  the  Sigma 
Meter  for  the  horizontal-path  laser  signal  and  the  differential  ther¬ 
mometer  were  digitized  at  15-s  intervals  by  means  of  a  Benson-Lehner 
LA.RR-M  digitizer.  The  resulting  data  were  written  on  magnetic  tape  to 
facilitate  later  use. 

F.  Results 

1.  Meteorological  and  Horizontal  Control  Path  Measurements 

Ground-based  measurements  were  obtained  on  three  days  during 
this  experiment,'  the  slant-path  measurements  were  limited  to  two  days 
because  of  the  problem  with  aircraft  availability.  The  results  of  the 
ground-based  measurements  are  illustrated  in  Figures  33-37.  Figure  33 
presents  scatter  diagrams  for  the  three  experimental  days  for  the  mea¬ 
sured  values  of  scintillation  (or  )  on  the  350-m  horizontal  control  range 

m 

versus  C  as  obtained  from  the  differential  thermometer  measurements, 
n 
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VERSUS  REFRACTIVE  INDEX  STI 
TURE  CONSTANT  (C*).  FOR  THE 
PATH  EXPERIMENT 
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FIGURE  34  TIME  VARIATIONS  ON  19  APRIL  1969  OF  (a)  LASER 

SCINTILLATION  ON  350-m  HORIZONTAL  CONTROL  RANGE 
<«-).  AND  tb)  THERMAL  TURBULENCE  LEVEL  {o(AT!| 
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FIGURE  35  TIME  VARIATIONS  ON  20  APRIL  1969  OF  (a)  LASER 

SCINTILLATION  ON  350-m  HORIZONTAL  CONTROL  RANGE 
<oml.  AND  lb)  THERMAL  TUPEULENCE  LEVEL  [o<AT)l. 
Downward  excursions  of  traces  were  caused  by  cloud  shadows 
drifting  across  range 
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The  analog  standard  deviations  usod  here,  as  obtained  from  the  Sigma 
Meter,  have  been  corrected  by  means  of  Eqs.  (13)  and  (15),  Section  III. 

The  dashed  curves  in  this  figure  represent  the  empirical  func¬ 
tion  fit  to  the  data  from  the  Woodland  horizontal  propagation  experiment. 
The  curve  fits  the  data  on  20  and  23  April  reasonably  well,  but  seems  to 
underestimate  the  measured  values  on  19  April.  The  reason  for  this  i** 
not  known,  but  it  could  be  due  to  an  error  in  calibration  of  tuc  amplifier 
for  the  differential  thermometer  on  that  day,  since  the  measured  values 

of  C  seem  abnormally  low  in  comparison  with  those  on  the  other  two  days, 
n 

when  the  general  weather  conditions  were  similar. 

A  number  of  the  data  points  on  the  three  days  are  anomalously 
high.  As  was  discussed  in  Section  III,  an  increase  in  scintillation  mag¬ 
nitude  occurs  when  the  laser  beam  drifts  and  the  receiver  is  no  longer 
located  in  the  center  of  the  beam.  This  is  what  occurred  at  various 
times  during  the  experiment.  When  this  condition  was  noticed,  it  was 
immediately  corrected  by  recentering  the  receiver  in  the  team.  The  effect 
of  this  realignment  procedure  is  apparent  in  Figure  34(a)  from  the  char¬ 
acter  of  the  trace  ct  about  1200  CST.  This  figure  illustrates  ihe  varia¬ 
tion  of  and  a(Ax)  with  time  on  19  April.  It  may  be  noted  that,  prob¬ 
ably  because  of  the  slightly  sloping  nature  of  thu  site,  the  transition 
from  unstable  to  adiabatic  to  stable  conditions  occurs  sooner  at  the 
location  of  the  differential  thermometer  (at  the  west  end  of  the  range) 
than  over  the  330-m  range  as  a  whole.  The  adiabatic  condition  is  indi¬ 
cated  by  the  minimum  In  d(«T)  at  1720  CST  and  the  d  minimum  at  1740  CST. 

m 

This  out-of-phase  condition  explains  the  reason  for  the  strange  shape  in 
the  envelope  of  data  points  near  the  origin  in  Figure  33(a). 

Figure  35  illustrates  the  time  variation  its  d  and  *(&T)  on 

m 

20  April.  Although  higher  frequencies  occur  in  the  latter  trace  because 
d( AT)  represents  a  point  measurement,  rather  than  a  line  integral 


does  a  ,  there  is  a  good  correlation  between  o  and  a(AT).  This  can  be 
m  m 

seen  more  clearly  in  Figure  36,  which  has  an  expanded  time  scale.  The 
dips  in  the  traces  were  caused  by  cloud  shadows  drifting  across  the  range 
from  time  to  time. 

The  time  variation  of  a  and  a (At)  on  23  April  is  presented  in 

n 

Figure  37.  As  indicated,  there  were  several  occasions  when  beam  drifting 
caused  the  receiver  to  be  off-center  and  necessitated  realignment.  No 
slant-path  data  were  collected  on  23  April. 

Figure  38  illustrates  the  vertical  temperature  profiles  obtained 
by  means  of  the  instrument  mounted  on  the  C-131  aircraft.  On  19  April 
at  1215  CST  there  were  temperature  inversions  present  at  approximately 
1200  m  (3900  ft)  and  2000  m  (6600  ft)  above  ground.  By  that  afternoon 
at  1520  CST,  the  two  inversions  had  apparently  merged  to  form  a  single 
inversion  layer  from  1700  to  2000  m  (5600  to  6600  ft).  Little  structure 
is  apparent  in  the  profiles  for  20  April,  except  for  isothermal  layers 
at  1500  m  (4900  ft)  at  0955  CST  and  at  1700  m  (5600  ft)  at  1320  CST. 

2.  Slant-Path  Moasurements 

For  a  slant  path  of  length  R  with  the  origin  of  coordinates  at 
the  receiver,  Tatarski  (1961)  has  shown  theoretically  that 

*  R  "11/2 

a  *  1.49  (2tt/\:7/12  f  C  (r)r5/6  .  (16) 

t  /  n 

-  o 

In  order  to  evaluate  the  integral,  it  is  necessary  to  specify  C  as  a 

n 

function  of  distance  along  the  slant  path.  For  our  work,  we  have  assumed 
that  C^  is  horizontally  uniform  and  decreases  exponentially  with  height 
(z)  according  to  the  model 
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C  =  C  e 
n  no 


kz 


(17) 


> 


where  k  is  a  selectable  parameter  with  units  of  m  and  C  is  the  near- 

no 

ground  (2  m)  value  of  C  .  This  function  is  graphed  in  Figure  39  for  two 

n 

reasonable  values  of  k. 


FIGURE  39  ASSUMED  EXPONENTIAL  FORM  FOR  VERTICAL 
PROFILE  OF  NORMALIZED  REFRACTIVE  INDEX 
STRUCTURE  CONSTANT 
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When  the  C  profile  is  substituted  into  Eq.  (16),  and  geomet- 
n 

rical  considerations,  involving  the  slant-path  elevation  angle  9,  are 
taken  into  account,  we  obtain  the  following: 

a  =  1.43(2tt/\)7/12C  A1/2  ,  (18) 

t  no 

where  the  range  integral  (a)  is  given  by 
Case  I  (ground-based  laser): 


R 

.  ,  -B(R-r)  5/6 

A  =  A.  =  /  e  r  dr 


(19) 


Case  II  (airborne  laser): 


R 

■/' 


,  -Br  5/6 

A  =  A  =  /  e  r  dr 


(20) 


and  B  =  k  sin  9.  The  value  of  the  constant  in  Eq.  (18)  has  been  changed 
slightly  from  Tatarski’s  plane-wave  value  to  incorporate  our  results  from 
the  horizontal  propagation  experiment. 

The  integrals  A  and  A  in  Equations  (19)  and  (20)  are  not 
1  2 

solvable  analytically  and  must  be  integrated  numerically.  Since  our 

experiment  involved  a  ground-based  laser,  Eq.  (20)  was  used,  and  values 

of  a  were  calculated  for  each  run  for  which  a  measured  value  of  o  on 
t  m 

the  slant  path  was  available,  using  values  of  C  derived  from  the 

'  no 

ground-based  measurements.  Values  of  o  were  obtained  from  a  total  of 

°  m 

28  runs;  the  remaining  10  runs  could  not  be  analyzed  because  oi  poor 
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3ignal-to-noise  ratios,  or  else  the  data  were  invalid  because  of  equip¬ 
ment  malfunctions. 

The  measured  values  (a  )  are  compared  with  the  theoretical 

m 

values  (c^)  in  Table  VI  and  Figure  40.  The  dashed  curves  in  the  figure 


Table  VI 


REDUCED  SLANT-PATH  DATA,  APRIL  1969  EXPERIMENT 


Run  No. 

0  (slant) 
m 

<°C) 
ff  (At) 

o (slant) 

lEBEKai 

1 

0.80 

m 

1.63 

1.17 

4 

0.94 

1.68 

1.17 

5 

0.87 

0.148 

1.63 

1.14 

6 

0.94 

0.156 

1.71 

1.20 

7 

0.71 

0.152 

1.44 

1.04 

8 

0.69 

0.140 

1.32 

0.95 

12 

1.10 

0.160 

1.73 

1.21 

13 

0.99 

0.120 

1.02 

0.71 

14 

1.13 

0.160 

1.72 

1.20 

15 

0.80 

0.120 

1.25 

0.90 

16 

0.64 

0.100 

1.04 

0.75 

17 

0.67 

0.108 

1.12 

0.81 

18 

0.62 

0.112 

1.04 

0.75 

19 

0.69 

0.100 

0.93 

0.68 

20 

0.62 

0.096 

1.01 

0.72 

21 

0.71 

0.088 

0.92 

0.66 

22 

0.92 

0.080 

0.87 

0.61 

23 

0.80 

0.076 

0.83 

0.58 

24 

0.85 

0.084 

1.00 

0.70 

26 

0.76 

0.080 

1.00 

0.70 

27 

0.71 

0.096 

0.95 

0.69 

28 

0.83 

0.120 

1.17 

0.85 

29 

1.29 

0.180 

2.39 

1.73 

30 

1.17 

0.200 

3.33 

2.39 

32 

0.87 

0.188 

1.82 

1.32 

36 

0.67 

0.180 

1.75 

1.27 

37 

0.67 

0.160 

2.10 

1.52 

38 

0.67 

0.048 

0.79 

0.57 

9  (THEORETICAL) 

TA-7472-41 

FIGURE  40  MEASURED  LOG-INTENSITY  STANDARD  DEVIATION  ON 

THE  SLANT  PATH  AS  A  FUNCTION  OF  THE  THEORETICAL 
VALUE.  ASSUMING  TWO  VALUES  OF  k  IN  THE  EXPONEN¬ 
TIAL  MODEL  Cns  Cn0  exp  (-k*).  The  dashed  curve*  represent 
the  empirical  function  previously  fit  to  the  data  from  the 
horizontal  propagation  experiment 


represent  the  same  empirical  function  previously  fit  to  the  data  from  the 
Woodland  horizontal  propagation  experiment.  Although  there  is  a  fairly 
large  scatter  in  the  data  points,  the  data  are  represented  reasonably 
well  by  the  empirical  function,  at  least  for  k  =  0.02  m  *. 

It  is  apparent  from  Figure  40  that  the  values  of  a  obtained 

t 

during  the  experiment  were  not  large  enough  to  provide  an  adequate  test 
of  the  supersaturation  effect  with  regard  to  scintillation  magnitude  that 
was  found  during  the  Woodland  experiment.  The  laser  we  used  for  the 
slant-path  experiment  did  not  have  enough  power  to  permit  measurements 
at  longer  ranges,  which  would  have  given  larger  values  of  o^.  There  is 
some  evidence  from  these  slant-path  data  that  scintillation  saturation 
is  occurring  in  the  range  from  ff  =  1,5  to  3.0,  but  measurements  at 
longer  ranges  will  be  necessary  to  confirm  the  supersaturation  effect  on 
the  slant  path. 

As  shown  in  Figure  40,  there  is  a  slight  tendency  for  the 
measured  scintillation  values  to  be  greater  when  the  propagation  path 
went  through  an  elevated  temperature  inversion  layer,  indicating  that  the 
thermal  turbulence  frequently  present  in  such  layers  may  be  significant 
in  enhancing  scintillation. 

Although  no  power  spectra  were  computed  for  t  slant-path  data, 
it  was  apparent  in  the  raw  data  (see  Figure  31,  Sectii  IV-E)  that  the 
scintillation  frequencies  were  much  higher  than  those  observed  during 
the  horizontal-path  experiment.  The  dominant  frequencies  of  scintillation 
on  the  slant  path  were  typically  in  the  range  1000-4000  Hz,  while  those 
in  the  Woodland  experiment  were  generally  on  the  order  of  100-200  Hz. 

These  higher  frequencies  are  caused  by  the  pseudo-wind  effects  caused  by 
the  movement  of  the  aircraft  and  laser  beam. 

Before  proceeding  to  the  next  section,  it  is  worth  noting  that 
Eqs.  (18),  (19),  and  (20)  predict  that,  for  the  same  range,  wavelength, 
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and  ground-level  value,  the  scintillation  magnitude  will  be  largest 
for  a  near-ground  horizontal  path  [for  which  Eq.  (18)  reduces  to  Eq.  (9)), 
next  largest  for  a  slant  path  with  a  ground-based  laser,  and  least  for 
the  same  slant  path  but  with  an  airborne  laser.  This  is  the  case  because 
thermal  turbulence  has  the  most  effect  near  the  laser  source.  The  equa¬ 
tions  also  show  that  for  a  given  range,  the  scintillation  magnitude 
decreases  with  increasing  slant-path  elevation  angle,  and  reaches  a 
minimum  for  a  vertical  (90°)  propagation  path. 
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V  INCORPORATION  OF  EXPERIMENTAL  RESULTS  INTO  IASER  EYE  SAFETY  GUIDELINES 


A.  General 

In  the  previous  report  on  this  project  (Johnson  et  al.,  1968),  a 
procedure  was  described  for  estimating  the  probability  of  eye  damage 
from  exposure  to  a  laser  beam.  Similar  work  has  been  reported  by  Deitz 
(1968,  1969).  The  new  eye-hazard  evaluation  procedure  presented  in  this 
section  incorporates  revisions  that  (1)  reflect  the  results  of  our  latest 
experimental  work,  and  (2)  streamline  the  procedure  to  make  it  easier  to 
use  on  an  operational  basis.  The  equations  and  assumptions  underlying 
each  part  of  the  eye-hazard  evaluation  procedure  (as  represented  by 
Figures  41-48  and  Tables  VII  and  VIII)  will  be  briefly  explained  in  a 
subsequent  section. 

It  should  be  emphasized  that  these  eye-safety  guidelines  should  be 
considered  tentative  and  used  with  caution  until  the  new  experimental 
results  that  they  reflect  can  be  confirmed  for  a  variety  of  laser-beam 
and  propagation-path  configurations. 

B.  Eye-Hazard  Evaluation  Procedure 

The  input  parameters  required  for  the  eye-hazard  evaluation  proce¬ 
dure  are  as  follows: 

<I>,  latitude  (deg) 

U,  surface  wind  speed  (knots) 

V,  atmospheric  visibility  (miles) 

N,  cloud  cover  (tenths) 

6,  slant  path  elevation  angle  (deg) 
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R,  range  <ka) 

P,  laser  peak  power  (W) 

Q,  laser  pulse  rate  (pulses/s) 

\,  laser  wavelength  (y) 

laser  bean  divergence  (mrad) 
d,  minimum  beam  diameter  (a) 

T,  exposure  tlae  (s) 

Laser  type  (CW,  long  pulsa,  or  Q-switched) 

Type  of  slant  path  (ground-based  or  airborne  laser) 

Date  and  tlae 

Using  the  above  information,  the  probability  of  eye  damage  any  be 
estimated  by  the  method  outlined  below. 

(1)  in  the  case  of  relatively  low-power  lasers,  go  to  Step  (2). 
For  other  lasers,  skip  to  Step  (3). 

(2)  Divide  the  LASER  POWER  P  (W)  by  the  minimum  beam  area 
Ttd2/4  (m2),  using  the  MINIMUM  BEAM  DIAMETER  d  (m)  to 
compute  the  maximum  power  density  available  for  eye 
damage  (W/a).  If  this  value  does  not  exceed  the 

appropriate  MAXIMUM  SAFE  POWER  DENSITY  LEVEL  (I  M  ) 

safe 

given  in  Table  VIII,  then  the  laser  can  be  considered 
safe  and  the  rest  of  the  procedure  disregarded. 

(3)  Use  TIME  OF  YEAR,  HOUR  OF  DAY,  and  UTITUDE  #  In 
Figure  41  to  find  SOLAR  ELEVATION  ANCLE  Or.  Interpolate 
between  graphs  for  dates  between  those  given.  For 
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latitudes  in  the  southern  hemisphere,  reverse  the  dates 
on  graphs  (a)  and  (c).  For  nighttime,  proceed  to 
Step  (5). 

(4)  Use  Of  and  CLOUD  COVER  N  (tenths)  in  Figure  42  to  find 
INSOLATION  CATEGORY  S  (strong,  moderate,  or  weak). 

(5)  For  daytime,  use  S  and  SURFACE  WIND  U  (knots)  in 

Table  VII  to  find  NEAR-GROUND  REFRACTIVE  INDEX 

-6  -1/3 

STRUCTURE  CONSTANT  C  (10  m  ).  For  nighttime,  use 

n 

last  two  columns  of  table,  given  N. 

(6)  For  a  near-ground  horizontal  propagation  path,  proceed 

to  Step  (8).  For  a  slant  path,  use  RANGE  R  (km)  and 

SLANT  PATH  ELEVATION  ANGLE  8  in  Figure  43  to  find  RANGE 
-ll/6% 

INTEGRAL  A  (m  ). 

(7)  For  a  slant  path,  use  A,  C  ,  and  LASER  WAVELENGTH 
\  (p)  in  Figure  44  to  fiud  LOG-INTENSITY  STANDARD 
DEVIATION  cr.  (Proceed  to  Step  (9)). 

(8)  For  a  near-ground  horizontal  path,  use  R,  C  ,  and  LASER 

n 

WAVELENGTH  \  (w)  in  Figure  45  to  find  LOG- INTENSITY 
STANDARD  DEVIATION  <?, 

(9)  Use  ATMOSPHERIC  VISIBILITY  V  (miles),  LASER  BEAM 

DIVERGENCE  8  (milliradians),  and  R  in  Figure  46  to 

-2 

find  NORMALIZED  MEAN  POWER  DENSITY  I/P  (»  )• 

(10)  Multiply  7/P  by  LASER  POWER  P  (W)  to  obtain  MEAN  POWER 

-  -2 

DENSITY  I  (W  *  ). 

(11)  U*o  \  and  LASER  TYPE  (Q-switched,  long  pulse,  uj*  CW) 
in  Table  VIII  to  obtain  MAXIMUM  SAFE  POWER  DENSITY 
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-2 

I  (W  m  ).  For  a  CW  laser,  use  the  following 

S3  ±6 

exposure  times:  (1)  <  1  ms  for  a  slant  path  in¬ 
volving  an  aircraft;  (2)  1-10  ms  for  static  paths 
and  U  >  10  knots;  (3)  10  ms  -  1  s  for  static  paths 
and  U  S  10  knots. 

(12)  Divide  I  „  by  1  to  obtain  SAFE-TO-MEAN-POWER- 
safe 

DENSITY  RATIO  I  ^  /I. 

safe 

(  13)  Use  a  and  I  /I  in  Figure  47  to  obtain  INSTANTANEOUS 
safe 

PROBABILITY  OF  EYE  DAMAGE  Y. 

(14)  Use  Y  and  EVENT  FREQUENCY  F  (s"1)  in  Figure  48  to  obtain 

INTEGRATED  PROBABILITY  OF  EYE  DAMAGE  T}  for  given  EXPOSURE 
TIME  T.  For  pulsed  lasers,  set  F  equal  to  the  LASER  PULSE 
RATE  Q  (pulses/s).  For  CW  lasers,  set  F  equal  to  the 
values  given  below: 


Static  path 


U  £  10 

knots: 

F  =  500 

U  >  10 

knots: 

F  =  1000 

Aircraft 

path: 

F  =  5000 
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FIGURE  42  INSOLATION  CATEGORY  IS)  AS  A  FUNCTION  OF 
SOLAR  ELEVATION  ANGLE  (?)  AND  TENTHS  OF 
OPAQUE  CLOUD  COVER  (N) 


Table  VII 


RELATION  OF  NEAR-GROUND  (2  m)  REFRACTIVE  INDEX  STRUCTURE 

CONSTANT  (C  )  TO  WEATHER  CONDITIONS 
n 


C  (10 
n 

-6  -1/3 

m  ) 

Nighttime  Conditions 

Surface 

Daytime  Insolation 

Thin  Overcast 

Wind 

Category  (S) 

or  »  5/10 

S  4/10 

(knots) 

Strong 

Moderate 

Weak 

Cloudiness 

Cloudiness 

<  3 

0.60 

0.30 

0.40 

0.30 

0.40 

3-6 

0.50 

0.40 

0.25 

0.20 

0.35 

6-10 

0.40 

0.25 

0.15 

0.10 

0.20 

10-12 

0.20 

0.12 

0.05 

0.05 

0.05 

>  12 

0.12 

0.05 

0.05 

0.05 

0.05 
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(X)  AND  A  NEAR-GROUND  HORIZONTAL 
PROPAGATION  PATH 


Table  VIII 


MAXIMUM  SAFE  POWER  DENSITY  LEVELS  (I  ) 

safe 


FOR  CORNEAL  EXPOSURES  (RECOMMENDED 
BY  USAF  MANUAL  AFM  161-8,  1969) 


I 


safe 


(W/m2) 


Laser 

Q-Switched 

Long  Pulse 

cw  I 

Wavelength  (a) 

Time 

(20  ns) 

(1  ms) 

<lms 

l-10ms 

10ms- Is 

1.0600 

Day 

3.21X10® 

0. 59x10° 

7.15X10^ 

... 

... 

... 

Night 

1.30X10* 

— 

— 

— 

0.6943  -  0.6328 

Day 

Night 

6 

5. 35XlOg 
0.98X106 

4 

1.43X10 

0.26X10* 

3 

2.86X10 3 
0.52X10 

1.43X10 3 
0.26X10 3 

7.15Xlof 
1.30X10 2 

0..145  -  0.4880 

Day 

— 

... 

3 

2.86X10 

1.43X10 3 

7. 15X10 2 

Night 

— 

0.52X10 3 

0.26X10 

1.30X10 2 

FIGURE  47  INSTANTANEOUS  PROBABILITY  OF  EYE  DAMAGE  (» 
AS  A  FUNCTION  OF  LOG-INTENSITY  STANDARD 
DEVIATION  to)  AND  SAFE-TO-MEAN  POWER  DENSITY 
RATIO 
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C.  Basis  for  Eye  Hazard  Evaluation  Procedure 

1.  Figure  41 

The  curves  in  this  figure  are  based  upon  the  following  equation 
(Smithsonian  Meteorological  Tables,  1963): 

<r  =  sin  1|sin®sinC  +  c0s®cos£cos[0.262(h-12)] |  ,  (21) 

where  a  =  solar  elevation  angle 

«  =  latitude 
£  =  solar  declination 
h  =  hour  of  day. 


2.  Figure  42 

\ 

\ 

The  relative  insolation  strength  (S),  on  a  scale  from  zero  to 
one,  is  given  by 


S  =  sin  a  (22) 

for  a  cloudless,  nonattenuating  atmosphere.  When  the  cloud  cover  (N)  in 
tenths  is  considered  and  the  albedo  of  the  clouds  is  taken  to  be  0.5, 

Eq.  (22)  becomes 

S  =  sin  [0(1  -  0.5N/10)]  ,  (23) 

which  is  the  basis  for  Figure  42.  The  categories  of  S  (strong,  moderate, 
and  weak  insolation)  were  determined  simply  by  dividing  the  range  of  S 
into  three  equal  intervals:  0  -  0.33,  0.33  -  0.67,  and  0.67  -  1.0. 

3,  Table  VII 

This  table  is  based  principally  upon  empirical  data  from  our 
experiments  and  those  of  other  investigators,  such  as  Wright  and  Schut* 
(1967),  Fried  et  al.  (1967),  Davis  (1966),  and  Goldstein  ct  al,  (1965). 
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Unfortunately,  most  of  the  values  of  C  reported  by  these  other  workers 

n 

have  not  been  accompanied  by  a  full  account  of  the  concurrent  meteorological 
conditions.  Hence  the  dependence  of  C  upon  weather  conditions  as  given 

It 

in  the  table  is  based  largely  upon  our  measurements.  The  values  in  the 
table  must  be  considered  only  as  approximate,  since  some  extrapolation  was 
required,  and  also  since  at  least  one  Important  variable,  ground  condition, 
has  been  omitted  for  the  sake  of  simplicity. 

The  basic  structure  of  the  table  is  patterned  after  that  given 

by  Pasqutll  and  Meade  (see  Slade,  1968)  for  the  estimation  of  atmospheric 

stability  categories.  Since  is  related  to  the  absolute  value  of  the 

stability,  appropriate  values  of  C  were  substituted  for  those  of 

n 

stability  category. 

4.  Figure  43 

Equations  (19)  and  (20)  (Section  IV-F)  were  used  for  the  compu¬ 
tations  for  the  curves  presented  in  this  figure.  The  value  of  the  constant 
k  was  taken  to  be  0.01  (see  Figure  39,  Section  IV-F)  in  order  that  the  eye- 
haxard  estimates  be  on  the  conservative  (safe)  side. 

5.  Figure  44 

This  figure  was  based  upon  Eq.  (18),  Section  IV-F,  for  the  slant- 
path  configuration,  but  we  have  limited  a  to  a  maximum  value  of  1.25.  The 
curves  are  thus  derived  from  the  Tatarski  relationship  but  take  into 
account  the  saturation  effect  indicated  by  experimental  results.  Since  the 
occurrence  of  scintillation  supersaturation  on  the  slant  path  has  not  yet 
been  adequately  confirmed,  we  have  not  made  any  allowance  for  this.  In  not 
doing  so,  we  have  erred,  if  at  all,  on  the  safe  side. 

6.  Figure  45 

This  figure  uses  Eqs.  (9)  and  (11),  Section  III-F,  and  incor¬ 
porates  the  supersaturation  effect  as  observed  during  the  horlxontal 
propagation  experiment.  It  is  significant  in  the  eye-hasard  context 
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that,  as  shown  in  the  figure,  the  scintillation  magnitude  (9)  reaches  a 
maximum  value  of  1.25  and  then  decreases  for  increasing  range  and  C^. 

7.  Figure  46 

This  figure  is  derived  piincipally  from  considerations  of  the 
laser-beam  geometry,  represented  by  the  equation 


4G 

i/p  =  -j-j 
tt3  r 


(24) 


where 


I  =  mean  power  density  (W/m  ) 

(assumed  to  be  uniform  across  the  beam) 

P  =  laser  output  peak  power  (W) 


0  =  beam  divergence  (milliradians) 
R  »  range  (km) 

G  a  atmospheric  transmission. 

The  transmission  is  given  by 


G 


-oR 

e 


(25) 


where  a  is  the  attenuation  coefficient  (km  *).  To  a  first  approximation 
the  relation  between  the  atmospheric  visibility  (V,  miles)  and  the  atten 
uation  coefficient  is  given  by 


V  =  2.43/a 


(26) 


for  visible  and  near-infrared  wavelengths. 


8.  Table  V1H 

.  2 

The  maximum  safe  power  density  levels  *,  m  )  given  in 

this  table  are  those  recommended  by  the  t'.S.  Air  Force  in  AFM  161-8 
(1989). 
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9.  Figure  47 


The  curves  in  this  figure  were  computed  by  means  of  the 
following  equations: 

go 

1  f  2 , 

V  =  /  exp(-x  /2)dx  ,  (27) 

M 

where  V  is  the  instantaneous  probability  of  eye  damage,  x  is  a  dummy 
variable,  and 

M  =  0/2  +  C^ln  (I  /I)  ,  (28) 

safe 

where  o  is  the  log  intensity  standard  deviation.  The  development  of 
these  equations,  as  detailed  in  the  last  report  (Johnson  et  al.,  1968), 
involves  the  assumption  of  a  log-normal  probability  distribution  of  in¬ 
tensity  fluctuations.  For  convenience,  these  derivations  are  reproduced 
in  the  Appendix. 

10.  Figure  48 

When  one  looks  at  a  laser  beam  continuously  for  T  seconds,  the 
integrated  probability  of  eye  damage  (r)  over  that  period  will  be  greater 
than  the  instantaneous  probability  of  eye  damage  (Y).  One  simple  way  to 
assess  this  increased  probability  is  to  assume  that  the  intensity  peaks 
occurring  at  a  single  observational  point  represent  separate  events. 

If  (l-Y)  corresponds  to  the  probability  of  no  damage  for  a  single  event 
(i.e.,  the  probability  that  any  one  intensity  peak  will  not  exceed 

I  ),  then  the  integrated  probability  of  no  damage  (l-F‘)  for  a  given 

sa  fr 

number  of  occurrences  (M)  in  the  time  period  (T)  is  given  simply  by 
multiplying  the  individual  probabilities: 
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i-r  =  (i-y)m  , 

or  T  =  1  -  (l-y)M 

Since  M  =  FT,  where  F  is  the  frequency  of  events, 

FT 

T  -  1  -  (l-Y)  •  (29) 

This  equation  forms  the  basis  for  Figure  48. 

The  event  frequency  F  may  be  considered  to  be  equivalent  to 
the  pulse  rate  for  a  pulsed  laser,  or  to  the  scintillation  frequency  in 
the  case  of  a  CW  laser.  If  we  assume  the  turbulent  elements  separated 
by  a  mean  correlation  distance  (Pq)  to  move  with  the  wind,  then  to  a 
first  approximation  we  have 

F  =  U  /p 
n  o 

where  U  is  the  wind  component  perpendicular  to  the  propagation  path, 
n 

Since  the  variation  in  p^  with  range  that  is  predicted  by  theory  dis¬ 
agrees  with  that  measured,  it  seems  best  for  our  purposes  to  assume  that 
F  depends  only  upon  wind  speed  and  beam  sweep  rate.  On  the  basis  of  our 
measurements,  the  following  values  have  been  selected  as  being  representa¬ 
tive  of  the  cases  given: 

Static  paths 

U  S  10  knots:  F  =  500 

U  >  10  knots:  F  =  1000 

Aircraft  slant  paths:  F  =  5000 

These  coefficients  allow  for  the  maximum  scintillation  frequencies 
observed  during  our  experiments. 
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11.  Setting  a  Level  of  Acceptable  Risk 


After  the  probability  of  eye  damage  (r>  is  found,  it  becomes 

necessary  to  compare  this  value  with  a  preestablished  value  (r^  ) 

representing  the  highest  probability  judged  to  be  an  acceptable  risk. 

The  magnitude  of  f  may  change  as  the  nature  of  the  situation  changes, 

-3  -5 

but  valuer  in  the  range  of  10  to  10  (1  chance  of  eye  damage  in  1000 

to  100,000  occurrences)  seem  reasonable.  By  way  of  comparison,  the 

probability  of  an  average  American  being  killed  in  an  automobile  accident 

-4 

in  the  year  1968  was  1  in  3600,  or  2.8  X  10  ,  while  the  probability  of 

-2 

being  injured  in  an  automobile  accident  was  1  in  100,  or  10  (World 
Almanac,  1970). 

D.  Examples  of  Use  of  Eye-Hazard  Evaluation  Procedure 

To  illustrate  the  use  of  the  procedure  given  above,  four  examples 
are  presented  in  Table  IX.  The  input  parameters  for  these  examples  were 
selected  to  be  as  representative  and  realistic  as  possible  and  to  cover 
a  variety  of  conditions.  Example  2  represents  the  conditions  of  an 
actual  recent  field  experiment. 
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Table  IX 


EXAMPLES  OF  USER  EVE-HAZARD  EVALUATION 


Parameter 

Example  1 

Example  2 

Example  3 

Example  4 

A.  SPECIFIED  PARAMETERS 

Date 

21  Dec 

1  May 

21  June 

21  June 

Time 

2200 

1200 

1500 

1500 

(LST) 

* 

30°N 

40°N 

20°S 

20°S 

N 

2 

0 

6 

6 

(tenths) 

U 

5 

8 

12 

12 

(knots) 

V 

20 

8 

20 

20 

(miles) 

Path  Type 

Slant 

Slant 

Horiz. 

Horiz. 

Laser  Site 

Airborne 

Ground 

— 

— 

0 

60° 

20° 

0° 

0° 

R 

0.4 

5 

1 

6 

(km) 

X 

0.4880 

0.6943 

1.06 

1.06 

(u> 

Q 

— 

0.5 

10 

10 

(pulses/s) 

0 

0.5 

0.5 

0.5 

0.5 

(mrad) 

d 

2 

- - 

— 

(lo’3m) 

P 

1 

7 

2.5X10 

7 

10 

7 

10 

(W) 

T 

1 

0.1 

0.1 

0.1 

(s) 

Laser  Type 

CW 

Q-sw. 

Q-sw. 

Q-sw. 

B.  DERIVED  PARAMETERS 

Step  No. 

3 

Of 

0° 

60° 

30° 

30° 

4 

S 

— 

Strong 

Mod. 

Mod. 

5 

c 

n 

A 

0.35 

0.40 

0.25 

0.25 

do'Vx/3] 

6 

6X103 

6X105 

— 

— 

(m11/6) 

7,8 

cr 

0.4 

1.25 

1.23 

0.80 

9 

i/p 

30 

-2 

5x10 

4 

-2 

9X10 

(m"2) 

10 

i 

30 

1 . 25x10® 

7 

4X10 

5 

9X10 

2 

(w/m  ) 

11 

Isafe 

1.3X10* 

5.35X10® 

3.21x10* 

3.21X10* 

(W/m2) 

12 

I  ,  n 

safe 

y 

4.3 

4.28 

8.0 

360 

13 

6X10-5 

•ixio"2 

io-2 

io'14 

14 

r 

-2 

2X10 

mo'1 

>io’1 

10- 12 

IS 

Safe  ?* 

no 

no 

no 

yea 

Based  upon  a  selected  value  of  T  »  10 


Appendix 

DERIVATION  OF  EQUATIONS  FOR  PROBABILITY 
OF  EYE  DAMAGE  FROM  A  SINGLE  EVENT 


Let 


where 


Now 


X  3  In  (A/A  ) 
o 


A  =  field  amplitude 

A  =  field  amplitude  reference  value 
o 


I 


where 


I  =  field  intensity  • 


Also,  from  Eq,  (A-l), 


hence 


A  =  A  eX  , 
o 


I 


2X 

e 


Now  we  assume  that  X  Is  a  normal  random  variable. 


to  assuming  that  A  A  follows  a  log-normal  distribution. 


shown  (Mood,  1950)  that  the  mean  value  of  I  is  given  by 


2  2X  i  2  /  -  2\ 

1  r  A  W  I  =  A  exp  I2X  +  2 

o  o  '  X' 


(A-l) 

(A-2) 


(A-3) 

This  is  equivalent 
Then  it  can  be 

(AH) 
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100 


as  given  by 
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ABSTRACT  (continued) 


either  a  ground-based  or  airborne  laser,  as  well  as  near-ground  horizontal 
propagation.  The  airborne-laser  slant-path  configuration  is  shown  to  give 
the  lowest  eye-damage  probability,  other  factors  being  the  same. 
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The  results  of  two  laser  scintillation  experiments  in  the  atmosphere  are  presented, 
as  is  a  revised  laser  eye-hazard  evaluation  procedure  incorporating  these  results.  TV 
objective  of  the  research  is  to  determine  the  influence  of  atmospheric  thermal  turbu¬ 
lence  on  laser  eye  safety.  » 

The  experimental  work  on  a  3.5-km  near-ground  horizontal  propagation  path  was  de¬ 
signed  to  determine  the  effect  of  range  upon  laser  scintillation  magnitude,  represent 
by  the  log  intensity  standard  deviation  (o).  It  was  found  that  for  very  short  ranges 

and  low  thermal  turbulence,  the  measured  values  (o  )  agree  well  with  theory  (a  ).  Hov 

ra  t 

ever,  for  longer  ranges  and  higher  thermal  turbulence,  the  values  of  a  reached  a  maxi 

mum  value  of  approximately  1.25  at  6^  =  2.50,  and  then  decreased  beyoncl,  fitting  the 

empirical  function  c  -  &  /(I  +  0.160  fairly  well.  This  "supersaturation  effect"  if 

m  »  t 

quite  significant  in  reducing  eye-damage  probabilities  for  typical  ranges  and  thermal 
turbulence  conditions.  j 

The  second  experiment  was  an  attempt  to  confirm  that  scintillation  saturation  occu 
on  a  slant  path  between  a  ground-based  laser  and  a  receiver  on  an  orbiting  aircraft. 
Some  evidence  of  saturation  was  found,  but  the  effect  was  not  clearly  established  be- 
I  cause  of  the  limited  ranges  necessitated  by  the  laser  used. 

The  eye-hazard  evaluation  procedure  is  quite  generally  applicable,  and  requires 
only  readily  available  input  data.  The  method  covers  propagation  on  a  slant  path  wit 
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